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CHAPTEB  I 

Introduction 

Recently  attention  has  been  directed  toward  the  effect  of  a 
magnetic  field  on  the  flow  of  an  electrically  conducting  fluid, 
usually  referred  to  as  magnetohydrodynamics  (MHD) .   Magnetohydro- 
dynamlc  flow  has  many  engineering  applications  such  as  raagneto- 
hydrodynamic  generators,  accelerators,  electromagnetic  flowmeter, 
electromagnetic  pumps,  and  similar  devices. 

Let  us  now  consider  that  the  electrically  conducting  fluid 

is  flowing  with  a  steady  state  velocity  V,  and  that  a  magnetic 

field  with  field  density  B   is  applied  perpendicularly  to  the 

flow  (Fig.  1) .   Because  of  the  interaction  of  the  flow  and  the 

magnetic  fields,  an  electric  field  E,  .   is  induced  perpendicular 

to  both  V  and  B  . 
-     -ap 

This  electric  field  is  given  by  the  following  equation: 

E,  ,   =  V  x  B   .  (1) 

— lnd     —   — ap  ' 

For  simplification  assume  that  the  electrical  conductivity, 
a   ,  is  constant  In  spite  of  the  magnetic  field.   By  Ohm's  Law  the 
current  density  induced  in  the  conducting  fluid,  and  denoted  by 


iind  i3i 


=iind  =  ae  Sind  •  <2> 


Simultaneously  occurring  with  the  Induced  current  is  the 
induced  Lorentz  force  F,  .  which  is  given  by  the  following: 

^ind  =  iind  x  2^  •  (3) 
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Fig.  I       Vector   diagram  of  MHD 


This  force  F,  .  occurs  because  the  conducting  fluid  cuts  the 
lines  of  the  magnetic  field.   Because  the  vector  product  of  equa- 
tion (3)  yields  a  vector  perpendicular  to  both  Jlnd  and  B   ,  the 
Induced  force  is  parallel  to  V  but  opposite  in  direction. 

For  the  mora  general  case,  we  further  consider  an  eleotrlo 

field  E   perpendicular  to  both  B   and  V,  but  opposite  in  dlrec- 
—ap  — ap     —       cr 

tlon  to  J«_d«  The  current  density  due  to  this  applied  electric 
field  is  «L,orld»  The  net  current  J  through  the  conducting  fluid 
is  then 

I     =   ae(Eap  ♦  V  x  Bap) 

The  ponderomotive  or  Lorentz  force  associated  with  this 
current  is  then 

F  =  J  *  S-„  =  <JJE=~  +  V  x  B  J  x  B_  .        (5) 
—     —   — ap      e  — ap   —   — ap    — ap  " ' 

In  equation  (5),  if  E   >  V  x  B   the  system  is  a  magneto- 
hydrodynamic  accelerator  (or  pump)  which  may  be  used  as  a  thrust- 
producing  device.   If  Eo„  <  V  x  B„  .  it  is  a  magnetohydrodynamic 
~ap    ~~   *~ap 

generator. 

The  equations  of  MHD  flow  of  continuous  fluid  media  are  the 
ordinary  electromagnetic  and  hydrodynamic  equations,  modified  to 
take  account  of  the  interaction  between  the  fluid  motion  and  the 
magnetic  field. 

On  the  assumptions  that:   (1)  the  flow  is  laminar,  (2)  all 
fluid  properties, f ,  0  ,  K,  ?  ,  are  constant,  (3)  the  displacement 
current  is  negligible  (as  in  most  electromagnetic  problems)  I.e. . 


no  oscillations  of  very  high  frequency  occur,  (H-)    the  permeabil- 
ity. yUet    and  conductivity,  a   ,  are  constant  scalar  quantities, 
and  (5)  the  effect  of  gravitational  force  is  negligible  and  the 
Lorentz  force  is  the  only  body  force  on  the  fluid,  the  basic  MHD 
equation  can  be  written  as  follows  (1,  2)+. 

Maxwell's  equations  in  cgs  electromagnetic  units: 

curl  H  =  J+nJ  ,         div  J  =  0  (6) 

curl  E  =  -yUe   ||  .     div  H  =   0  .  {?) 

Ohm's  law  for  a  moving  fluid: 

J  =  >„e(E  +  V  x/^H)  .  (8) 

Continuity  equations 

div  V  =  0 

The  modified  Navler-Stokes  equation: 

^   +  (V  grad)V  =  -J   grad  p  +  j>v\   +  A  (J  +^eH)    (9) 

For  steady-two-diraenslonal  flow  with  the  usual  Prandtl 
boundary-layer  assumptions  that: 

s2u  a2u 

(1)   — p  is  very  small  compared  with  — £  so  it  can  be 

s*  <?y2 

neglected; 


*  Numbers  in  parentheses  refer  to  references  at  the  end  of  the 
chapter. 


(2)  the  transverse  velocity  v  is  small  in  comparison 
with  uj  and 

(3)  consequently  the  pressure  gradient  §§  is  a  function 

^x 

of  x  alone, 

equation  (9)  is  resolved  into  the  following  three  equations  in 
the  x,  y,  and  z-dlrections,  respectively: 

ay2 


»S"8  -   -i|f^y*f  U/z-W  •  <10, 


o  =   -SS*^.wa-W  ■  (11) 

0   =  ^e(JxHy  -  JyHx)  .  (12) 
Next,  Ohm's  law  may  be  written  as: 

Jx     =     °e(Ex  */W    •  I") 

Jy     =     °e(Ey  "^V    •  U*> 

Jz     =     °e(Ez  ^e^y  WV    ■  <x5) 

If  the  magnetic  field  is  steady,  j#  =  0,  and  from  equation  (7), 
curl  E  =  0,  that  is 

•B 

(16) 

(17) 

(18) 
Equation  (6)  becomes: 
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k„Jz 
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5A 

ax 

ay 

(19) 
(20) 

(21) 

For  the  simplification  of  the  foregoing  equations  further  assump- 
tions are  made  as  follows:   (1)  Variations  in  the  z-direotion  are 
assumed  to  be  zero,  (2)  the  electric  field  term,  E  ,  measured 
across  the  insulated  duct  walls  Is  zero,  but  small  local  values 
may  exist  in  the  midstream  region;  however,  these  will  be  con- 
sidered negligible  and  E  is  taken  as  zero  at  all  points  across 
the  duct.   This  implies  J   Is  also  zero. 

Thus  taking  Ey  =  0,  J  =  0,  and  j^  =  0,  In  Ohm's  law  and 
Maxwell's  equations,  it  can  be  shown  that 

E   =  constant  =  EQ  , 

B   =  constant  magnetic  field  applied  =  BQ  , 

H2  =  Ex  =  Jx  =   0  . 

Equations  (10),  (11),  and  (12)  reduce  to 

■S**8  -     -/fi^-^O  +  uBo-^)  .      (22) 


If     =     °eBx(E0  +  ^0   "  vBx>    •  (23) 

0     =     0.  (2*0 


Equation   (21)    may  be  written  as, 

WEo  *  ^o  "  VV  =  "  FT  •  (25) 

Unlike  the  case  of  the  fully  developed  velocity  profile,  the 
foregoing  equations  are  coupled  by  the  appearance  of  v,  the 
transverse-velocity  component.   Since  the  applied  external  mag- 
netic field  has  no  x-  or  z-component,  and  B  =  0,  it  follows  that 

the  component  of  the  magnetic  induction  in  the  x-directlon,  B  , 

x 

is  induced  only  by  J  .   It  is  assumed  that  B  is  negligible  in 
comparison  with  the  applied  field  BQ.   These  assumptions  reduce 
the  number  of  the  equations  to  two: 

5  u   av 

^  +  7?     =  °  '  (26) 

u£u  au  _  1  dj>  afu  _  fgfa  , 

ax  T  v  ay      f  dx  +  v       2    /»   Vilo  *  ^o'       lZ7) 

which  become  the  basic  governing  equations  for  magnetohydro- 
dynamlc  channel  flow. 

In  this  report  we  consider  a  laminar  flow  of  a  conducting 
fluid  with  constant  properties  entering  a  semi-infinite,  non- 
conducting flat  duct  with  a  normal  transverse  applied  magnetic 
field.   The  initial  velocity  profile  is  assumed  to  be  uniform  or 
a  nonmagnetically  fully  developed  parabolic  profile,  and  at  a 
large  distance  downstream  from  the  entrance,  the  velocity  profile 
is  the  fully  developed  Hartmann  profile.   There  is  also  an  exter- 
nal variable  resistance  connecting  the  two  perfectly  conducting 
end  plates  (electrodes)  which  are  displaced  to  infinity  (Fig.  2). 


o 

CO 


CD 


CO 


GO 


00 


CD 


■*- 

I 

o 

o 

E 

■4*1 

CD 

c 

y 

y 
y 

y 

y 

c 

(/)        . 

y 

y 

a 

>           ° 

y 

y 
y 

-C 

y 
y 

y 

y 

y 
y 

o 

c      > 
O 

y 

y 
y 

_D 

-J=    E 

y 

y 

Q_ 

i_ 

^ 

y 
y 

-C      o 

y 

y 

y 
y 

"ai 

■»  1 

y 

y 

7~~ 

y 
y 

y 
y 

o 

j     / 

y 

i_ 

c      T3 

y 

y 

y 
y 

I 

P\ 

y 
y 

y 

o 
a. 

.2       C 

y 

>~ 

y 
y 

o 

CD 

y 

y 

y 

a> 

y 

y 

>N 

<y     •— 

y 

y 

/             '             \ 

y 
y 
y 

1_ 

c 

y 

i                                                \ 

y 

CD 

E 

o     o 

y 

y 

i_      — 

y 

y 
y 

i 

> 

m 

V 

1     1 

1      1 

, 

y 

y 

o 
a; 

OJ      c 

e 

3 

CO 

en 

In  steady,  laminar,  non-magnetohydrodynamlc  flow  In  the  en- 
trance of  the  channel,  boundary  layers  form  at  both  walls  due  to 
viscous  friction.  Near  the  entrance,  these  boundary  layers  form 
much  like  those  for  flow  over  a  flat  plate.   However,  in  channel 
flow  the  restriction  of  the  flow  near  the  walls  is  compensated  by 
an  Increase  in  the  rate  of  flow  in  the  center  of  the  channel. 
The  boundary  layers  then  Increase  in  size  as  the  flow  progresses 
down  the  channel,  until  they  merge  into  each  other  to  asymptotic- 
ally form  the  Poiseuille  velocity  profile. 

In  the  case  of  flow  of  a  conducting  fluid  under  the  effect 
of  a  transverse  magnetic  field,  i.e. .  the  magnetohydrodynamlc 
flow,  the  boundary  layers  will  form  on  both  channel  walls  in  a 
similar  manner,  but  the  force  due  to  the  magnetic  field  adds  an 
effect  to  the  flow.   As  in  the  non-magnetic  case,  the  thickness 
of  the  boundary  layers  will  gradually  increase  as  the  velocity 
profile  asymptotically  approaches  the  Hartmann  fully  developed 
profile. 

In  many  magnetohydrodynamlc  applications,  e.g..  MHD  genera- 
tors, liquid  metal  pumps,  plasma  engines,  etc.,  the  flow  of  fluid 
is  seldom  fully  developed  and  is  of  the  boundary  layer  type  with 
variable  longitudinal  pressure  gradients.   For  this  reason  a 
study  of  the  velocity  fields,  boundary  layer  development,  and 
friction  factors  in  the  entrance  region  of  MHD  channel  is  of 
practical  Importance  and  has  been  a  subject  of  investigation  in 
recent  years. 

In  this  entrance  region  we  wish  to  obtain  expressions  for: 
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(1)  the  pressure  drop  between  any  two  sections,  (2)  the  velocity 
distribution  at  any  section,  and  (3)  the  value  of  x  at  which  the 
fully  developed  flow  is  attained. 

The  two  simultaneous  equations,  equations  (26)  and  (27)  with 
appropriate  boundary  conditions,  are  sufficient  to  solve  for  the 
unknowns,  u  and  v.   The  pressure  ceases  to  be  an  unknown  function 
since  it  can  now  be  evaluated  from  the  potential  flow  in  the 
central  core  by  the  aid  of  the  Bernoulli  equation. 

Although  equations  (26)  and  (27)  are  much  simpler  than  the 
original  equations  representing  governing  equations,  their  exact 
analytical  solutions  have  not  been  found  so  far.   However,  a 
number  of  different  approximate  solutions  has  been  published  in 
channel  entrance  flow  problem.   Hwang  and  Fan  have  published  a 
bibliography  of  hydrodynamlc  entrance  region  flow,  a  thorough  re- 
view and  classification  of  the  literature  on  the  subject  which 
has  appeared  in  the  last  hundred  years  (3). 

The  approximate  solutions  may  be  classified  In  four  general 
categories  --  the  momentum  integral  method,  linearlEatlon  method, 
matching  method,  and  finite  difference  method. 

In  the  momentum  integral  method  the  flow  is  divided  into  a 
boundary  layer  part  near  the  wall  and  a  potential  flow  part  In 
the  central  core.   A  parabolic  velocity  profile  (or  any  other 
similar  velocity  profile)  is  assumed  in  the  boundary  layer  and  is 
joined  with  the  center  core  velocity  profile  which  is  assumed  to 
be  a  straight  line.   A  momentum  integral  equation  is  derived 
based  on  the  momentum  conservation  principle.   This  approach  was 
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devised  and  applied  by  Schiller  C+)  for  flow  in  a  circular  tube 
and  is  similar  to  the  Karman-Pohlhausen  momentum  integral  method 
applied  to  a  flat  plate.   Schiller's  method  was  applied  to  the 
magnetohydrodynamic  flow  by  Maciulaitis  and  Loeffler  (5).  by 
Moffatt  (6)  in  1964,  and  by  Tan  (7)  in  1965. 

In  the  second  category  of  solutions,  the  inertia  terms  of 
boundary  layer  equations  are  linearized.   This  category  of  solu- 
tions is  capable  of  providing  continuous  solutions  for  the  veloc- 
ity of  distribution  and  pressure  drop  in  the  entrance  region  of 
nonmagnetic  flow.   This  class  of  solutions  for  circular  tubes  is 
mainly  due  to  the  work  of  Langhaar  (8,  9)  In  19^0  and  Targ  (10) 
in  1951.   This  category  of  solutions  was  applied  to  magnetohydro- 
dynamic flow  in  the  entrance  region  by  Hsueh  (11)  in  1963  and  by 
Snyder  (12)  in  1965. 

The  third  group  of  solutions  is  constructed  by  matching  the 
boundary  layer  solutions  which  are  valid  near  the  entrance  with 
the  perturbations  of  the  fully  developed  solutions  which  are 
valid  for  downstream.   This  class  of  solution  was  originally 
given  by  Schlichtlng  (13)  in  193^  for  a  flat  duct.   This  matching 
method  was  applied  to  magnetohydrodynamic  flow  at  the  entrance 
region  by  Barsness  (14)  in  i960  and  by  Roldt  and  Cess  (15)  in 
1962. 

The  fourth  approach  involves  reduction  of  the  continuity  and 
momentum  equations  to  finite  difference  equations  which  are 
solved  numerically  on  an  electronic  digital  computer.  This 
method  was  used  by  Bodoia  and  Osterle  (16)  in  1961.   Hwang  and 
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Fan  (1?)  in  1963,  Shohat  et  al.  (18)  in  1962,  and  Hwang.  Li  and 
Pan  (19)  in  1966  applied  this  method  to  MHD  channel  flow. 

The  purpose  of  this  report  is  to  study  laminar  magnetohydro- 
dynamlo  flow  in  the  entrance  region  of  a  flat  duct  by  presenting 
in  detail  the  solutions  according  to  the  momentum  Integral  method 
(5),  Targ's  linearization  method  (12),  and  the  matching  method 
(15). 
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NOMENCLATURE 

a  Channel  half-height 

B  Magnetic  induction 

e  Electric  field  magnitude  factor,  EQ  ■  -  eu0B0 

E  Electric  field  intensity 

H  Magnetic  field  intensity 

J  Electric  current  density 
M 

p  Fluid  pressure 

u  x-component  of  velocity 

v  y-component  of  velocity 

x,y,z  Space  coordinate 

P  Kinematic  viscosity 

f  Fluid  density 

M  Dynamic  viscosity 

jj  Magnetic  permeability 

a  Electric  conductivity 


Hartmann  number,  M  =  Ba(-^r)* 


Ik 
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CHAPTER  II 


Matching  Method  (1) 

In  this  analysis  we  consider  the  laminar  flow  of  an  electri- 
cally conductive  fluid  with  constant  properties  entering  a  semi- 
infinite,  nonconducting  flat  duct  with  a  normal  transverse 
applied  magnetic  field.   The  fluid  is  assumed  to  enter  the  duct 
with  a  uniform  velocity  profile.   There  is  an  external  variable 
resistance  connecting  the  two  perfectly  conducting  end  plates 
(electrodes)  which  are  displaced  to  infinity. 

The  process  of  analyzing  the  velocity  field  in  the  entrance 
region  is  patterned  after  a  method  developed  by  Schllchtlng  (2) 
for  non-MHD  flow  in  which  the  flow  Is  divided  into  two  sections: 
flow  near  the  inlet  or  upstream  section,  and  a  downstream  section 
which  approaches  the  fully  developed  flow.   In  the  upstream  sec- 
tion, the  flow  is  analogous  to  boundary  layer  flow  over  a  flat 
plate  with  a  pressure  gradient;  hence,  the  solution  is  modeled 
after  the  Blaslus  series  solution.   In  the  downstream  section, 
the  velocity  distribution  is  assumed  to  be  the  sum  of  the  fully 
developed  profile  and  a  deviation  of  the  profile  from  its  asymp- 
totic parabolic  distribution.   The  integration  is  then  performed 
in  the  upstream  direction  in  order  to  find  this  deviation  veloc- 
ity.  After  obtaining  both  solutions  in  the  form  of  series  ex- 
pansions, they  are  Joined  at  the  point  where  both  solutions  are 
valid.   In  this  way  an  approximate  description  of  the  flow  field 
in  the  entire  magnetohydrodynamic  entrance  region  is  obtained. 

For  the  magnetohydrodynamic  flow  the  governing  equations  are 
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C 

o 

CD 

o 


c 

CD 


CD 

£ 

^.    c 
<->  £ 

-o  +- 
o 

S  E 


c 

D 


a; 


2  = 
!=;   ° 


c    +- 

o  -r 

U    5 
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(see  Chapter  I) : 

ax   *  ay  "  u  ■  u> 

and 

3  y 

Upstream  Solution 

The  upstream  section  Is  divided  Into  the  subreglons  of  ln- 
vlsold  potential  flow  and  viscous  boundary  layer  flow.   In  the 
potential  flow  where  u=U,  v=0,  and  U  Is  only  a  function  of  x, 
equation  (2)  becomes 

u  dx  -    f   dx    /   U     f  l3J 

Solution  of  equation  (3)  for  g*  an(i  substitution  into  equation 
(2)  yields 

9  y 

Let  the  transformation  variable  J  be 

?   • 

and  assume  the  potential  flow  velocity,  U,  to  be  in  the  follow- 
ing form: 

U(x)   =   u0(KQ  +  Kx|  +  K252  +  K3^3  +  ...)  (5) 
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In  which  KQ  is  equal  to  one,  since  U  =  uQ  at  ?  =  0,  and  K, ,  K-. 
...  are  unknown  constants  to  be  determined. 

Let  the  other  transformation  variable  f,   be 


and  assume  that  the  stream  function,  /'(x.y),  can  be  expanded  Into 
a  series  In  )j  In  a  form  similar  to  that  for  U(x) .   Thus, 

where  s  is  the  width,  and  the  f  are  functions  of  >7  only.   To 
satisfy  continuity  equation,  equation  (1),  u  and  v  are  given  by 


37   '  v  =   " ax 


(7) 


From  equation  (5)  we  have 

U  51  =  u0(1  +  Kl?  +  K2^  +  K3*3  *  •••>  (  U0(K1  *  2*K2 
♦  352K3  +  ...)}■ 


'£ 


where 


df 

1            P 

-      1         ^x      /x 
2           2         ' 
V  ^UQ 

i  y 

dx     " 

'      2,2 

V  s   uQx 

-      2   x    * 

Therefore,  we  have 

U  dx"  "  -2T  (1  ♦  Kl?  ♦  K2^  +  •••)(K1  +  2JK2  +  352K3  ♦  ...) 

u2r 

-     "^  (K!   +    (K2  +  2K2>5  +   3(KXK2  +  K3)$2  +    .....J.  (8) 
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From  equations    (6)    and    (?)    we   have 

u(x.y)      =    §1    =     u0s(|f.(t)    +  &&)   +   ...)f*. 
where 


3K  1        ^0  1        f_^0     ,  _1_ 

?y     =      2V>»x     =      2  ,/  p  x      /s     =      2sj  * 

Therefore,    we  have 

u(x,y)      =     ^  (f'W)   tff'lf)   +  J2f}(?)    +   ...  j 
Similarly,    from  equations    (6)    and   (?)    we   obtain 


(9) 


=     -uQs  [txW   *  2Jf2(t)   +  3?2f3(r)   *  -..} 


ii 

dx 


-u0s{jq(7)  ♦  J2f$<?)  +y>?yi)  ♦  ...}  y 


where 


dS        i  5 

dx     "     2  x 


ax_2Vy2x  J_  2x    * 

Therefore,    the   transverse  velocity  v(x,y)    becomes 

un5s 


v(x,y)      =     -  -1^-  (f^f)    ♦  2   f2(?)   +   352f3(?)   ♦   ...] 
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From  equation  (9)  we  obtain 


X 


uo5 


(f£(>r)  +  2^f3(r)  +  352f^(?)  f  ...j 


T^-{fj(f)  *rf|(f)  +  f 2f'3(?)  ♦ ...) .  (id 


It  =  iKfi<*>  +5f2^  *i2q^)  *  ».}/? 

=   — ■£-%{  f^'(f)  *  jff^'(r)  ♦  52f3"(r)  +  ...]   .         (13) 

Substituting  equations  (5)  through  (13)  into  equation  (^), 
we  obtain 

2 
^  (qtjr)*  f^(f)+|2f3(r)  +  ...){^(?)*2Jf3(?)+352fi(?)-»-...] 

2 

■  "Hf(  fi(?)+Jf2^)','52f3^)-,-"-){fi(?)',-ff2(?>^2f3^)*-"} 
2 

■  sf  {fi(r>+  ^(r)+|2f3(n*...){f1(r)>2|f2(t)+352f3(l?')+...j 
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^u0 
~8x~ 


^  2 


-J^  {  %  ♦  (K2  +  2K2)  f  f  3(^2  +  K3)|  2  «•  . .  A 

2 

♦  ^-^{(l+Klfi.K252+K353+...)-  |[q(?)  +  Jf2(Z) 

+  ,)iffi,^*/f2,(?»+A3,(?)t-) 


8scf 
where 


2  2  2 

a  BAu.  a  B.a    , 

e  0  °  _  (  e  0   wJi)  u 

f          ~  ffi          a2   ° 

-  «2  <4,)<^)  uQ 


2   "J2u0 


and 


Du0 

- 

VUQ 

2 

u0 

8s2J2 

8s2    ^       ' 

3    UQ 

-   ox 

Substituting  these  into  the  above  equations  and  simplifying,  we 
obtain 
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-{f»(f)+5f«(?)+52f»(^)+...}(f1(^)+2jf2(^)+3|2f3(r)+...} 

=  kf^   +  (K2  +  2K2)f  ♦  3(^2  +  K3)52  *  ...j 

+32M2J2{(l+K;J+K252+K3f3+...)-  |(fi(t)+Jf2(f)+52f^(f )♦...)] 

♦  {fi'(?)+Jf2,(?)+52f3,(?)  +  ---  }  •  (1*0 

Equating  like  powers  of  ?  gives 


-0 
for  5  , 


fj'(?)  f  f£(?)f(?)   =  0  (15) 


for  f1. 


for  5 


f{(t)^(1T)  -  2q(f)f2(f)  -  f^(?)f1(^)  -  ry(?)     =  4KX  (16) 

2 


f22(?)+2q(7)f$(<o  -  3f3(r)q(ir)  -  2f2(f)fg<r)  -  ^(fc^ 

=      8K2  ♦   i+K2  +  f»'(?)    ♦   16M2(2   -   fj(?))  (17) 


for  if3. 


3q(?)f^(r)  ♦  3f2(?)f5(?)  -  wj(t)f^dr)  -  3f3(f)f^(?) 

-  2f,,(?)f»(?)  -  fx(f  )fjt(tf) 
=   12K3  +  12KXK2  +  f£'(f)  ♦  16M2(2K  -  fgif))  (18) 

where  M  =  (ae^pj?)*  BQ  |  is  the  Hartmarm  number.   The  boundary 
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conditions  for  equations  (15)  through  (18)  are  determined  from 
u=0  and  v=0  at  the  wall  (y=0),  and  become 

fx(0)   =  f2(0)   =   ...   =  fn(0)   =   0 

f£(0)   =   f£(0)   =   ...   .  f£(0)   .   0 

At  the  edge  of  the  boundary  layer,  % t*>  °°   and  u  =  U,  so  that 
equating  equations  (5)  and  (9)  gives 

u0(KQ  +  Kxf  +  K252  +  K3f3  ♦  ...) 

=  -£■  f{(t)  +jf2(?)  +  fzm?)  ♦  ... 
t{(7)  *$t,z(i)  +f2r'3(?)  +j3f;(r)  t  ... 

=   2  +  2K1f  +  2K2^2  ♦  2K353  ♦  ...  (19) 

Collecting  coefficients  of  like  powers  of  f    yields 
f{(»)   =  2 

f2(»)   =  2KX 


fA<»>   =  2Kn-l  • 

For  large  •?,  the  foregoing  relations  may  be  Integrated  to  yield 
the  approximate  relations 
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f2  ~  2K^  -f  A2 

*  •  • 

•  •  • 

fn  ~  «ft-ll  *  An  ' 
Further,  from  continuity,  one  obtains 

a 
uQa  =   /   udy  =  [^(x,y)] 

0  y=a 

At  the  center  line  of  the  duct,  1    is  large  and  equations  (20) 
are  valid.   In  addition,  at  y=a  one  finds  that 

*.  J-   ^ 


(20) 


uos2 


Using  these  relations  in  the  foregoing  equation  together  with 
equations  (6)  and  (20)  yields 

f  (x.y)      =      2uQa  (jffj^)  ■+  ^2f2(f)    ♦  %  *tjft)    *    . . .) 

1  =  2&h  +  Ai»  +|2(27*  V  *j?3<zf*  V  *  •••} 


or 


or 


1     =      1  +   2JA1  ♦  KXJ   +   2A25<:  +  ^K2  +   2A353  ♦    ...    , 
and   thus, 
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Kx     =   -2A1  ,     K2  =   -2A2 Kn  =   -2AR  . 

Substituting  these  into  equations  (20),  we  obtain 

K 


f 


n  =  2Kn-l?--F 


for  large  values  of ^  .   The  K  -values  needed  in  each  equation  can 
be  determined  f ronf  the  solution  of  the  preceding  equation  of  the 
set  (15)  to  (18),  and  the  entire  upstream  flow  field  is  thus 
described. 

Downstream  Solution 

After  Schlichting,  the  velocity  in  the  downstream  region  Is 
taken  as 

u  =  uf  (y')  +  u'U.y')  ,  (22) 

In  which  uf(y')  and  u' (x,y')  are  the  fully  developed  velocity  and 
deviation  velocity,  respectively. 

In  the  region  of  fully  developed  flow,  the  momentum  equation 
may  be  written 

0  -    1  dp     *2u   °eB0  u   °eB0E0 

f  dx +  y 37z        r  r  (23) 


and  the  solution  by  Hartmann  (3)  is 

_  dp/dx  +  qeBQE0      cosh(My,/a) 

f         "  R2      KL  coshM   ' 

°eB0 

°e  i 
where  M  =  («y)   BQa,  the  Hartmann  number. 


(24) 
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Integrating  equation  (24)  from  zero  to  a  with  respect  to  y' 
and  remembering  the  assumptions  that  dp/dx,  a   .  BQ,  and  EQ  are 
constant,  we  obtain 

a  dp/dx  +  oBnE.   a  r     cosh(^J-), 

n2  Ly  McoshM   slmU    a   JJ  „ 

aeB0  ° 

a(df  *  cBoV    TM   -   tanhMl  ,,,., 

^2  I  «  J'  <25) 

The   continuity   integral 

a 

/     u-.dy«      =     u.a 

0        *  U 

can  be  Introduced  into  the  left  hand  side  of  equation  (25)  to 
reduce  it  to  the  form 

a(df  *   aeW  f  H  -  tanhMl 


u0a  =        „2 1  "  -M- 

°eB0 


from  which  we  can  obtain 


5J  ♦  °eB0E0        U0M 


Q2     _  M  -  tanhM  * 
°eB0 

Substituting  the  last  expression  back  into  equation  (24)  we 
finally  obtain 
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u.M   r     cosh(M  *-) I 
Vy,)   "  M  -  tanhM [  1  '    coshM J"  (26) 

When  M  is  equal  to  zero,  this  can  be  written  as 

u.M  r  cosh(M  X—)  i 

Llm     uf(y)      =       Llm     M   .   tanhM    [l   -        co3hM J 

n  'm  *U  1*1  ■    *u 

u.M[ooshM   -   oosh(M  ^-)] 

=        Lim     . 

M K>         [M   -   tanhM  ][coshM] 

This   can  be   evaluated  according  to  L'Hospltal's   rule   as   follows 

u-M [coshM   -   cosh(M  ^-)1 

Llm     -2-= ^-L- 

M— *0         [M   -   tanhM]  coshM 

u0[ooshM-cosh(M  2-)]    +   u0M[slnhM-(^-)sinh(M  ^-)] 
_        Lim     _ 

M K)  (1-sech  M)coshM   -    (M-tanhM) sinhM 


Lim 


u0[sinhM-(il^-)2sinh(M  ^-)]  +u0[sinhM-(^-)  sinh(M  £■)]    + 
M *0        2sechMsechMtanhMcoshM  +  (1-sech  M)slnhM  +  (1-sech  M) 

uQM[coshM-(^-)2cosh(M  ^-)] 
sinhM  +  (M-tanhM) coshM 

un[coshM-(^-)2cosh(M  £-)]  +un[coshM-(^-)2cosh(M  ^)]  + 

L,im  r r, P 

M K)  2sechJM-2tanhMsechMtanhM+2sech  MtanhMslnhM+( 1-sech  M) 

un[coshM-(J!^-)2cosh(M  ^1)1  +U.M  [sinhM-(^-)  3sinh(M  ^-)1 

U  L  ft S   J U   ^ 8.  &  ■* 

coshM*2sech2M2tanhMslnhM+(l-sech2M)coshM+(l-sech2M) 

coshM+(M-tanhM) sinhM 


Lim 
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=    —  - 

That   is 

uf(y)      =     |u0[l-<£)2]  (27) 

which  is  the  equation  for  the  velocity  profile  for  non-magnetic, 
fully-developed  channel  flow. 

For  the  deviation  velocity  u'(x,y').  Schllchtlng  assumed  a 
series  of  the  form 

4\n12 


u'(x.y')   =  u0Cie   1   <p'(y'/a) 


-*Wi  s2 
*   u0C2e   25  dj£(y'/a)  +  ...  (28) 

where  the  constants  C  and  k,    the  f  unction  <f(y'  /a) ,  and  the  region 
over  which  (28)  holds  are  to  be  determined.   In  the  downstream 
section  the  deviation  velocity  is  never  very  large  and  hence  the 
higher  order  terms  in  (28)  may  be  neglected  as  an  approximation. 
Thus, 

u'(x.y')   =  u^e    L      (£>'(y'/a)  .  (29) 

Substituting  equation  (22)  into  (2),  we  have 


30 


,  lW9tt'i  ,duf       sa1, 

(u^U'H^)    ♦    V^   +  jyr) 


.2  _ 

f  dx+    ^v"    ,2  +  "    ,2J 

J  dy'         ay' 


°eB0uf        aB0U*        °eB0E0 


(30) 


f  ?  t 

Differentiating  equation  (30)  with  respect  to  y' ,  the  pressure 
term  vanishes  because  it  is  a  function  of  x  alone. 

daf  an'   aul^ul   ,   .  n  a2u' 
dy'  ax  *  ffp  dt    *  luf+u  '  axay1 

wd3uf    a3u->    °eBo  duf    3U-,    !a!aifa„n 

where  — — j-  =  0  since  EQ  is  assumed  to  be  a  function  of  z-direction 
only. 

From  equation  (22),  we  have 

S  -  £■  •  <32> 

And  according  to  the  continuity  equation,  ^—  ■►  —7  =  0,  equation 
(32)  becomes 

—  -   -  —  (33) 

py«   _    ax  '  Kii> 

Substituting  equation  (33)  into  equation  (31)  and  simplify- 


lng,  we  obtain 


U-fU') 


32u-      ,dS-   92u., 
dy*    sy' 


Prom  equation  (29),  we  have 


ay' 


3X  ay'  ax 

?~  e 

ay'-  a' 

.2 


.3_uJ.  Vl      ""M 

,.3     -       J 


From  equation   (26),    we   have 


duf  uQM2  slnh(^) 
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'&-£S» -=£*♦»!• 


-lA  's2 
-UllJf2  -^,$2 

=    -IT"  W  Pi<V>   •  (35) 

-TT  e  Ji<V  (36) 

*2u«  4VlGl^2      -^l^2  V' 

-  ax  e  *1^>  (37) 


-,2    ,  unC,      -4L5' 

3   u'  0  1  ,,       l'     r  „,/y*>  /-Q. 

~72     =     —T  e  £l   'a5  (38> 


^  3d  e   ~1'     ylYjXL,    .  (39) 

<?y  *  a     a 


dy'      =      "   a(M-tanhM)    *      ooshM      '  ^°^ 
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d2u„  unM3  oosh(^-) 

dy'2     "  a2(M-tanhM)  C0ShM      ' 

d\  uV+  sinh(^-) 

7^     =     "  a3(M-tanhM)    (      ooshM      }    "  ^ 

From  equation    (33)  ■    we   get 
V     =     -    Jf^dy. 

-4-unX.,Cn52      -iJ-X.,?2  . 

4Vici?2   -^i^2  ,      v' 

4au  X  C  52     -la,?2 

=        %        e  9i^)  •  <w> 

And  from  the  transformation  variable 


Vx 


V  *\ 


we  get 


1?=  —   •  (14) 

Substituting  equations  (28)  and  (3D.  and  equations  (38) 
through  (4-5)  into  equation  (36)  and  neglecting  the  higher  order 
terms,  we  obtain,  > 
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U   -   tanhM   (1   "      coshM      )    +   u0Cle  <Pi^>J 

*[-^.«-^W>] 

*  x 6  *1(T>  L  "   a2(M   _   tanhM)    <    cosh  R      ' 


+   - 
a 


2      e  Tl    *  a  ;J 


=  ~^_i'a2(M-ta„hM)    (    coshM      '  a3  ^(a)J 

"      f      [~   a(M   -   tanhW)    l    oosh  M      '  a  T  lv  a  '  J 

In  which,    in  the  last   term, 


4.* 

a 

«2      ^2a2"c 


2 
a 


H52«2 


3^ 
Thus,  the  previous  equation  can  be  rewritten  as  follows 

ax     e      'Pi1  a  J\«-tanhH  L     cosh  M   J  +  ule     YVa'j 
Ifu -X ^C,!2   -^i?2,   v,  /   M3     cosh(^)       -)u?2    v,  \ 

*    *ll    9    *  Pi<V»(ii=fegr  coshM    -  v      *>!'<*•>} 

ax      IM   -   tanhM  I     oosh  M      J        ul  il    'a'] 

i*u2S2M2  ,  2  f  slnh(S£-)  ,  -^x52  „,    i 

*~^E—  (-M-   tanhwl      coshM      J  *   Cle        '      *!<¥>}•  (45» 

2   2 
Dividing  both   sides   of    the   above   equation  by  ^u„5  /ax  an<i   sim- 
plifying  yield 

HV  i^tfnhM  I  I"      cosh  I     fe'^lCle  $  S^J  'f  <¥> 


2  M  v1  ? 

+   ,   c   e"^   -JSL-f^ftLu  (IL,-,   A"8^'*    (^)-f)-(^) 
*   ^l8  M-tanhML   cosh  M      Jfl1  a  '    Xl°le  Tll  a  '  *1    la' 


.2 


In  the  above  equation,  second  order  terras  of  (p  can  be  neglect- 
ed.  This  is  the  same  approximation  made  by  Schlichtlng  (2)  in 
the  solution  of  the  non-magnetic  problem.   Also,  dividing  both 
sides  of  the  above  equation  by  \-,C-,e       we  finally  obtain 
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M\,    r    coshf^) 


il  *  \.  M  -  tanhK  L     cosh  M   J     /  *  1 

M  -  tanhM  I  cosh  »   Ju  v   ' 


with  X-i  as  Its  eigenvalue.   The  boundary  conditions 


v  =   0 

and 

gu 

■ay 

= 

0 

at 

y 

= 

v  =   0 

and 

u 

= 

0 

at 

y 

= 

are  equivalent  to 

$1  =  ££  =  0     at     y  =  0               C*8) 

^Pl   =  ^l   =   °      at      y   =  ±  a  .              (49  ) 
The  other  boundary  condition  is  chosen  as 

$>{  a  1     at     y  =  0  .  (50  ) 

This  is  possible  because  C.  is  still  free.   Having  these  boundary 
conditions,  equation  (4-7)  can  be  solved  by  a  power  series  method. 

Upstream  Pressure  Distributions 

Near  the  entrance  of  the  duct,  the  boundary  layers  are  very 
thin.   Most  of  the  flow  is  composed  of  the  accelerating  potential 
flow  in  the  central  core  of  the  duct.   Therefore,  the  pressure 
distribution  can  be  obtained  from  equation  (3), 
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_  1  d£      dU  "A  ^ofo  (  . 

f  dx  -   u  dx     f        f   '  *51' 

In  the  above  equation,  according  to  the  definitions 
^2a2uQ 


and 


Eo     = 

"eu0B0    ' 

one   can  obtain 

< 

„2   2 

aflQa 

>P 

f       " 

r^   " 

2 
a 

= 

M2  4 
a2 

i|M2?2u0 

and 


f 

= 

aBoeuo 

f 

= 

2    2 

oBoa        uo 

f  *      '    a2 

= 

2   *   S 
-   VT   -    . 
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Also  by  definition  we  have 
dx   _   2  x 


dx   =  f   df  . 

Hence,  according  to  equation  (5).  equation  (51)  can  be  rewritten 
as: 

-  j   dp  =   UdU  +  8u2M2f  (1+K^2  +  K2^3  +  Kj$k   +  .  ..)dj 

-  8^u2M2?d5  . 

Integrating  the  last  equation  from  x  =  x.  at  which  P  =  P.  and 
f =  0  to  any  section  in  the  upstream,  we  obtain 

p  «  e 

-  i  /        dp     =     i  /      d(U2)    +   8u2M2   /        (1+Kn52  +   K,53 
•>    PQ  ^   0  u        0  12 

♦   K  ^  +    ...)df  -   8eu2M2   /$f  d£ 

J  u        0 


P      -  P 
°f        =    |  "ouo  *  Ki^  ♦  K2^2  *  K3?3  ♦  •••)2 


?  ?     ?2  ^3  W*  e5 

3u0M      (  2    ♦  Kl  T  ♦  K2  T  *  K3     5    +    ' ' ' } 


3u2eM2    (■!-)    . 
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Finally,  the  pressure  drop  may  be  expressed  as: 

P   -  P 

jS _  m     z^  +    (2K2  +  K2  +  8M2  [l-e])?2 


J-  f    • 
2  ?U0 


♦  (2K3  +  2K1K2  -t-  y  A^3 

+  (2K^  +  2KXK3  ♦  K2  +  ^M2!^)^  ♦  ...     (52) 


Downstream  Pressure  Distribution 

At  the  end  of  the  upstream  section,  the  potential-flow  re- 
gion is  quite  small  and  the  acceleration  of  the  remaining  core 
flow  in  the  downstream  section  is  assumed  to  be  negligible.   This 
permits  the  use  of  the  momentum  equation  for  the  region  of  fully 
developed  flow 

o  _  .ifc+  „_!« .  !i!fc  .  !!ofo 

t   dx      y2     F  S 

in  the  downstream  section.   By  taking  a  force  summation  over  a 
differential  volume  of  fluid  according  to  the  above  equation  we 
can  obtain 

2 
1  dp.       £   ,^L\  aeB0u0    °eB0E0 

"  f   dx  -  "a  l ay' 'wall  *   f        +   f 


dp       TV   .  au  .         „2       D2 
dl  =   "  T  ^wall  *  °eB0u0  "  °eBOeu0 
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-fi    =     -T^W'^oV1-)    ■  '53) 

This   is  due   to   the   fact   that   the   pressure   varies   with  x   only  and 
■rj  can   be   evaluated   at   the   walls,    y'    =  +   a.      According  to   the 
solution  for  the   downstream  velocity  profiles   equation   (53)    can 
also   be  written  as   follows 


12/  uQM2  slnh(^)  ] 

a    \"  a(M   -tanhM)  I     cosh  M      J 


Vi    -^52 


1       ?l<£>wall}dX 


^M2f2u2 

*   (l-e)dx 


#,„  m2   sinhM  liv  e2 


2 

a 


0  f  coshM      ^   „    ■         l'     ti,„,,12i    ,f 

~  I  "  H   -    tanhM  ♦   Cle  iPl(1)JTdf 


iM2?2u2e 

j-2-   .    &   (l-e)dj 

^uQ        Wu^2,       w2   tanhM  -**!**  , 

"^5 fp I  -   M-tanhM      *   Cle  H<«)*f 

UM2u^?(l-e)5  d$ 


IK) 


2 

-  dp     =      8fu0   M   .   tanhM  gdf  *   8C1fu0  $»(l)e  £    df 

+   8M2/u2    (l-e)^df  (54) 

Letting  the  value  of  £  at  the  Jointing  point  be  denoted  by  ^., 
and  integrating  equation  (54)  from  ?=  ?t  at  which  P  =  P, ,  to 
|;  =  1y  at  which  P  =  P,  one  can  obtain 

£  -^     =     8^0  M^nT  '/*    *$  *   BClfu^»(l)    /«  e"^2r  d5 
1  1  *1 

♦   8M2fu2(l-e)    /*£d5 


Pl   "  P 


\?< 


(M_tanhM_    v    M       ,,.2l,,2        -2. 
U    -    tanhM   ♦    a"e)M     }V       "*l> 

2C,f?(l)        -Wf2  -"HUT? 

*  — ^j* (e        1        -   e        *   X) 
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as  the  pressure  distribution  in  the  downstream  section. 

For  the  complete  solution  of  the  pressure  distribution,  only 
C1  and  3^  remain  to  be  evaluated.   These  values  are  determined  in 
the  process  of  matching  the  upstream  and  the  downstream  solu- 
tions. 

Matching  of  the  Upstream  and  the  Downstream  Solutions 

The  upstream  and  the  downstream  solutions  of  velocity  and 
pressure  distribution  are  to  be  matched  at  a  point  where  both 


solutions  are  valid.   Because  there  are  two  undetermined  quanti- 
ties, C,  and  1;, ,  in  the  equations  to  be  joined,  two  conditions 
are  needed  to  solve  for  thera.   One  condition  is  based  on  the  fact 
that  the  centerline  velocities  represented  by  both  solutions  are 
equal  at  the  Joining  point.   The  other  condition  Is  that  the 
slopes  of  the  pressure  distribution,  j^,  for  the  two  solutions 
match  at  the  Joining  point.   These  two  conditions  are  used  to 
develop  two  simultaneous  equations  from  which  C.  and  *f ,  can  .be 
determined.   It  is  a  reasonable  assumption  that  the  upstream 
solution  is  more  accurate  than  the  other j  hence,  this  solution  is 
extended  as  far  downstream  as  possible  before  Joining.   Kj,  is  the 
highest  order  coefficient  solved  for  in  the  potential  flow 
velocity  expansion  and  it  is  assumed  that  when  KJ?   is  equal  to 
5  per  cent  of  U/UQ  in  equation  (5)  the  ^  term  can  be  safely 
neglected.   Thus 

U  =  u0(l  +  K-^   +  K2?2  +  K3£3  +  K^)  .  (56) 

The  downstream  centerline  velocity,  from  equations  (28)  and 

(3D  is 


u 


uf (y*)  +  u'(x,  y1) 


^n?2 


uf(0)  *  uQCe    1      (J)J(O)  . 


According  to  equation  (12)  and  the  boundary  condition 
represented  by  equation  (50),  the  above  equation  can  be  written 
as 


k2 


U  =   M  -  tanhM  (1  "  c7ihM>  *  Ve        '  ^l) 

Equating  the  centerllne  velocities  represented  by  equations 
(56)  and  (57)  gives 

I  ♦  Kjf  ♦  K2f2  +  K^3  +  K^4 

=   M  -  tanhM  (1  "  3oThM)  *  Cle  (58) 

The  upstream  pressure  slope  equation  is 

If  =   "  ^   g  "  °eB0U  "  °eB0E0  ' 

Substituting  the  centerllne  velocity  into  the  upstream  section 
and  using  the  definitions 

2  2 

2     aeB0a 
M   =   '  «,-. 

?i> 

and 

E0   =  '  "  eu0B0  ' 
we  can  obtain 

g  =   -fu2(l  t  K1?  +  K252  ♦  K^3  ♦  K^)(K1  ♦  2K2? 

t-  3K3f 2  t  «^3)  ff  +  M2  ^^  (1  ♦  Kx5 

a 

♦  K2§2  +  Kjf 3  *  K^)  ♦  M2  ^^  . 

Since  ^   =  -*-s —  and  -r—  =  pT,  the  upstream  pressure  slope 
^a  uq 


ita 


equation  can  be  reduced  to 


+   2K25  +   3K3?2  ♦   ^53)    -  M2(l  ♦  Kxf   +  K252 

+  K^3  +  K^)    +  M2e  ].  (59) 

On  the  other  hand,  from  equation  (53) t  the  downstream  pressure 
drop  is 

dp     fS>   ,8U,         Q2   ,,   , 
S  =  T  'FnU  "  aeB0u0(1-e) 

tl   i"*   »U»,       geB0u2   f^u0  ,-   , 
"   a  W  +  ay' 'wall  "   j»P   *  ~^5   l±"e' 

=  ^{"i^tin^^v"^^)-^-)}      Ceo, 

Equating  equations  (59)  and  (60),  we  have 

jfe  (1  +  Kx?  ♦  K252  +  K3?3  +  K/+5'4')(K1  +  2K25  +  Ji.*? 

+  *HC^3)  -  M2(l  +  K^  ♦  K252  -t-  K3?3  +  K^) 

2  -Ux  e^ 


M   tanhM 
M  -  tanhM 


-^n?2 


+  M2  -  C1e   1   <p£(0)   =   0  (6L) 


Eliminating  C^e       from  equations  (58)  and  (61)  gives  a 
seventh  order  equation  in  $ 
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(1  +  K.J  +  K2?2  +  K3f3  +   K4r)(Ka  +   2K^  +   3K3?2  +   4K^3) 
-  8M2(?  ♦  K,f 2  +  K2^3  +  v*  +  ^5)   +  8^™  5 

+   8M^  +   B   -TanhM   "   -   cl^h¥>  *I(1>?   -   8(1  +  Kl*   ♦  K2^ 
+  K^3  ♦  K^)£fJ(l)5     =      0 


v    /  ow    ^2   8M2  tanhM  g  ,     8M     .,     1  u„m 
H  *{  ^Z   *  "-1  *  M  -  tanhW  ?  +  M  -  tanhM  u  "  coshl-l'  $1VJJ 

-  8f!"(l)  5  +   3K3  +  31^  +  K1  -  8M2K1  -  8K2f^(l))52 

+  ['4-K1  +  ^K1K3  +  2K2  -  8M2K2  -  8K2<^(l))f3 

+  (SK^  +  5K2K3  -  8WK4  -  8K3  ^  ( 1 )}  5  ^ 

+  (6^^  +  3K2  -  8M2K^  -  8K^<j,£(l)}f5  ♦  {7K3Kj56 

+  (<+K2}^7  =  0  (62) 

Solving  this  equation  for  the  various  values  of  M  gives  f  at  the 
Joining  points.   Substituting  the  values  of  5  into  the  equation 
for  the  upstream  oenterline  velocity  gives  the  values  of  C.  for 
each  magnetic  field  strength  calculated. 
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NOMENCLATURE 

a  Channel  half-height,  a>| 

B  Magnetic  field  intensity 

e  Electric  field  magnitude  factor,  E.  =  -  euB. 

E  Electric  field  strength 

fn  Blasius  coefficient,  function  of  ^ 

J  Electric  current  density 

M  Hartmann  number,  M  =  BQ  a(^j)* 

p  Fluid  pressure 

a  Channel  height,  s  =  2a 

u  x-component  of  velocity 

uQ  Fluid  velocity  at  inlet 

uf  Fully  developed  velocity 

U  Potential  flow  velocity 

v  y-component  of  velocity 

w  z-component  of  velocity 

y  Direction  component  measured  from  channel  wall 

y'  Direction  component  measured  from  channel  center  line 

1  Transformation  variable,  1  =  i    ,  — 

I       2  J  i>x 

k  Eigenvalue 

P  Kinematic  viscosity 

?  Fluid  density 

a  Electrical  conductivity 

Y  Blasius  stream  function 

(p  Deviation  velocity  function 

f*„  Wall  shear  stress 
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^     Transformation  variable,  f=  (    ?) 

U0S 
Subscripts 

0  =    constant  magnitude 

1  =    Joining  point 

2  =    entrance  length 
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^7 
CHAPTER  III 

Schiller's  Method 
In  the  analysis  by  momentum  Integration  (1),  we  also  con- 
sider a  laminar  flow  of  a  conducting  fluid  with  constant  proper- 
ties entering  a  semi-infinite,  nonconducting  flat  duct  with  a 
normal  transverse  applied  magnetic  field.   In  this  case,  the 
velocity  profile  is  approximated  by  a  curve  comprised  of  two 
parabolas  and  a  straight  line  (see  Fig.  1);  the  vertices  of  the 
parabolas  lie  on  the  border  of  the  boundary  layer.   The  parabolic 
velocity  profile  in  the  boundary  layer  is  assumed  to  be 

§  =  !-[<*>  "I]2-  (1) 

where  U  is  the  velocity  of  the  freestream  and  is  a  function  of  x, 
the  coordinate  along  the  direction  of  the  flow;  6  is  the  thick- 
ness of  the  boundary  layer. 

In  this  study  two  different  boundary  conditions  are  con- 
sidered. One  is  the  uniform  velocity  distribution  across  the 
entry  plane,  and  the  other,  the  velocity  at  the  entry,  is  the 
fully  developed  non-MHD  one  which  has  a  parabolic  profile.  Both 
develop  in  the  entrance  region  of  the  channel  to  the  Hartmann 
velocity  profile  which  is  the  fully  developed  velocity  of  MHD 
flow. 

First  of  all,  the  momentum  integral  equation  for  flow  in  the 
duct  is  derived.   Consider  the  flow  region  within  the  boundary 
layer  between  section  ab  and  cd,  which  are  dx  apart  (see  Fig.  2). 
Let  the  x-axls  be  in  the  same  direction  as  the  flow,  and  let  the 
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y-axis  be  measured  from  the  wall.   Then  the  mass  flow  rate  Into 
the  section  ab  within  the  boundary  layer  is 

6 
Min  =  J"   2u  x  2w  x  /  dy 

6 
=  >+w  f  S     udy  . 

0 

The  mass  flow  rate  out  of  the  region  is 

6  H    6 

Mout  =  4w  f   ;   *+*»fst/   udyJ dx  ' 

The  difference  between  M,   and  M   .  is 

in      out 


H 


d 


6 


...  -Vt  =  -^stj  udy]d* 


which  is  the  mass  flow  rate  into  the  region  from  the  center  core. 
Therefore,  the  momentum  flux  into  the  region  due  to  the  bulk 
flow  (or  convective  flow)  is 

-tufO^f/  udy]  dx  +  4w  f  /  u^dy  . 
ax   0  0 

The  momentum  flux  out  of  the  region  due  to  the  bulk  flow  is 

*«f/     u2dy  +»wfi[J     u2dyl  dx   . 
0  dx     o 

Thus   the  net  momentum  flux  out  of   the   region  is 

4w  f  dx"  t  ;     u2dyjdx  tlwf  U^[/S   udyjdx   .  (2) 
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For  steady-flow,  the  net  momentum  flux  out  of  the  region  Is  equal 
to  the  net  force  acting  on   the  region  in  the  same  direction.   The 
forces  exerted  on  the  region  are  the  pressure  force,  the  shear 
force  of  the  wall,  and  the  Lorentz  force.   Hence,  when  Hall  ef- 
feots  and  induoad  magnetic  fields  are  assumed  to  be  small  enough 
to  be  negligible,  we  get 

SFx  =   "  2rw(2w)dx  +  l*r(/   JBdy)dx 

+  P(26)2w  +  (P  +  £   |2  dx)  .  2w  .  2d6 

"  (P  +  df  dxH2a  +  2d6)2w  (3) 

where  0  <  £  <  1.   In  equation  (3),  the  first  term  of  the  right- 
hand  side  represents  the  shear  force  of  the  wall.   The  second 
term  is  the  Lorentz  force  in  which  j  is  the  current  density;  B  is 
the  imposed  magnetic  field  intensity.   The  third  term  is  the 
force  acting  on  the  boundary  layer  at  section  ab,  the  fourth 
term  is  the  force  acting  on  the  boundary  layer  from  the  free- 
stream,  and  the  last  term  is  the  force  acting  on  the  boundary 
layer  at  section  cd. 

Rearranging  and  neglecting  the  infinitesimal  terms  of  higher 
orders  yields 

6 

SFx     =      -  /j-wr    dx  +   l*w(/      JBdy)dx 

0 

-  *»*  a£  dx  •  eo 

According  to  the  momentum  theorem,  equation  (2)  Is  set 
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equal  to  the  right-hand  side  of  equation  {k)    and  we  oan  obtain 
the  balance  of  forces  and  momentum  fluxes: 

ttffcll     u2dyjdx  +   tafU^[/6   udyjdx 

=     4»r,  dx  +  4w(/      JBdy)dx  -  4sw  £2  dx 
0  ax 


V«  If  -  J  Js*y  =   -  f «  £  5  ud*  -  f&  J  ^ 


(5) 


The  freestream  momentum  equation  is 

6  ""  "fi^J  (6) 

and  Ohm's  law  is 

J   =  <j(E  -  uB)  .  (7) 

These  provide  the  means  for  expressing  the  pressure  gradient  and 
the  Lorentz  force  in  terms  of  velocities  and  the  magnetic  field. 
That  is,  a  rearrangement  of  equation  (5)  and  the  substitution  of 
equations  (6)  and  (7)  into  equation  (5)  give 

?w  =   "  f  U  &  Ady  -  fA.   /6u2dy  +  ;6JBdy  .  6  |£ 


f U  £  /   (U-u)dy  -  f U  A  /  Udy 


6 

r 

0 


+  -Pd^/   (uU-u)dy  -  fA  /  uUdy 
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+   /      <J(2-uB)Bdy 
0 

dU 


s[-fug+  o(B-UB)B] 


-      CU^  /6    (U-u)dy   -   6  fug 

tf|rJ     (uU-u2)dy  -  f  4-0  /     udy 
ax   Q  ax        0 

+  OEB   /      dy  -   aB26U   /     7T  d(£> 
0  0      U        5 

*  6  U  fx  "  6EB<J  *  6°Ub2 

=  f  0  &  /   (U-u)dy  ♦  j3  i  /   (uU-u2)dy 
ax  Q  ax  Q 

+  60UB2  -  50UB2  /  $  d(£)  . 
0   U    6 

Introducing  the  displacement  thickness  6*  and  the  momentum  thick- 
ness 9  defined  by 


6*U      = 

6 
/      (U*-u)dy 
0 

and 

2 
9U        = 

6 
/      u(U-u)dy   , 
0 

respectively,    reduces   the  equation  into  the   form 

7W     =      &oBZU 

fl-/X§d(f>]  ♦•fujg^f; 

5^ 


In  order  to  solve  this  equation  for  the  boundary  layer 
thickness,  we  use  the  assumed  velocity  profile  expressed  in 
equation  (1).   That  is,  it  is  assumed  that  the  velocity  profile 
consists  of  the  parabolic  variation 

across  the  boundary  layer  and  a  uniform  core  velocity. 

First  of  all,  we  can  substitute  ^(— )   ....  for  the  shear, 
which  is  evaluated  from  the  assumed  velocity  profile.   Thus 

1  w  ~  ^l«'. 


'ay' wall 


y=0 


2MM 


ff-]y= 


y=0 


=  f^U  .  (9) 

The  terms  on  the  right-hand  side  of  equation  (8)  can  also  be 
evaluated  as  follows; 


/X&d(Xj   =  ?[1   -{(f)  -l)2]d(*) 


[(f)2  -*<f>3]" 


3  V  Jc 
3  (10) 
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;6  -  (i  -  -) 


=  J6[^)  -  5(f)2  +  Mf)*  -  <?>V 

"  l  6    3  62    63   57J0 


=   6 


[l.f+l-i] 


dx 


iLi  _   2a  (4i  (11) 

15  "   15  V  l 


2a  A  /A)  n2) 

15  dx  V  vi^ 


The  core  velocity  Is  allowed  to  vary  In  the  x-dlrectlon  in 
order  to  satisfy  the  conservation  of  mass  equation  applied  across 
the  entire  channel  cross  section.   If  u  is  the  average  velocity, 
the  principle  of  conservation  of  mass  gives 

6 
'J-awu  =  l+w  /  udy  +  4-wU(a-6) 
0 

=  ^wlJ  /  [2(£)  -  (£)2]dy  +  4wU(a-ft) 
0     5     6 

=   »HrtJ  [  i-  -  i  **]   +  i+wU(a-6) 
6    ">  6   0 

=  ^wU6(|)  +  ^wU(a-6) 
=   ifwUa(l  -  i  |) 


Thus, 


dU     _     UJL  1 

dx  dx   (1   _  1  A) 
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-    -     i  -  4  &  nil 

U  3a-  ll-,; 


3  a' 

(£)  (it) 


3(1    -if)2    dX       a 

f£(U29>     =     fu2ff.2fueg. 

From  equations   (11)    through   (It)    this   can  be  valuated  as 

f-*-   (U29)      =     f  li2    .    f|  f   (4) 
dx  '  15  dx      a 

+   2fU   •    —   (^)    ^ —  (£■) 

*   *ru        15    V       ,  1  4>2   dx   V 

-ni   "   3   a' 
9oPii2    ^  4aufU(£)  , 

=     J6[l-2(f)    +   (f)2]dy 

■     [6  -  6  +  f]    =     3   •  <l6> 
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Substituting  equations  (9).  (10).  (12),  (U),  (15).  and 
(16)  Into  equation  (8)  we  obtain 


H  -   *B2»a[l-f].|fU3(i  _\^2A(|) 


15   dx  V  +  45(1  _  1  k)Z   dx  laJ  * 

By  multiplying  both   sides  with  ^j  and  noting  that   (6^/)*Ba  =  M, 
kauf//4  =   Re,    and  u/(l   -  i  A)    =  U,    the   above   equation  becomes 


=>2lt.2 


^aufU(^)2 


3  a' 


-  d(£) 
90(1  -  ±  £)Z/"U       a 


aB   a  V       .  Jffau 

-   — 2 —  dx    ■ 


(*>2 


M 


^ 


-^  d(4)   +  *Ml 

72(1  -li)2     V         - 
3   a 


ttt  d(a)  +  ^-  • 


60(1  -li) 


a(*)2 
90(1  -  i  f)2       a' 


,A,2  M2  ,  Bea<|)2         „A.  aea(|)  , 

(o)      -7-  dx  +  '1    .    ?   d(£     +  § — a— _  d(£) 

72(1   -  i  A)2        a  60(1   -  £  i)        a 

3  a' 


T 


3  a' 


H   a(i) 
e    xa' 


4.2 


90(l-if)2d(a) 


=     (i) 


2  M* 


a-      T1* 


Bea  675(|)2+810(&)(1-  J  A)+5»Q(&)2 


18600(1   -ii)2 
3  a' 


<) 
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=   (!)2^dx+R 


7(f)2  +  6(f) 
360(1  -  i  A)2 


d(|) 


dx 

aH_ 


7(A)2  ♦  6(A) 


•2  6_2 

2 
a 


36o(i  -V^Xl  -U)2 


(17) 


Equation  (1?)  Is  to  be  integrated  from  x=0  to  x  with  the  boundary 
conditions  for  the  uniform  entry  velocity,  when  x=0,  6=0  and  when 
x=x,  6=6.   Thus 


a 


=*-  /  ^  =  / 


7(^)2  +  6(A) 


aH 


e  0 


36ofl-(^)(A_)][l-i(f)] 


2d<a"> 


(18) 


In  order  to  carry  out  the  integration  on  the  right-hand  side  of 
equation  (18),  an  integration  by  partial  fractions  is  performed. 

Let 

7(A)2  ♦  6(A) 


360[l-^(A)2j[i.i(A)]2 


"5*  va'    Jl  *        3xa' 

A(A)    +   B 


h{ 


36"o 


[i-^(f)2J      f1"!^      5-1 


l  -i(A)W 


35o  I-  : „2    :  ~ : Tz — 


[i-¥cMfi.j<i)r 


°fr  -  V  d>2J 


}• 
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Then  we  have 


A(f)  -  §  A(£)2  ♦  |  A(f)3  +  B  -  §  B(f)  +  i  B(£)2  ♦  C 

-  J  0(ft)  -  Jf  0(A)*  ♦  ^  .  i  0(A)3  +  d  -  ^  ^ 

=  7(f)2  ♦  6(A)  . 

Equating  coefficients  of  like  powers  of  (•*■)  in  the  above  equa- 
tion we  obtain 


B  ♦  C  +  D 


0    . 


A   -  |b   -   j  C      =      6, 


2a+1b_MIg 

3A*9BT> 


IA  +  £C    _     0 


2 
M 


7   . 


Solving  these  simultaneous  equations,  we  obtain 
3i\-)   [lM  ♦  162(^-)] 

A   =  g   *      , 

2 

63  +  891(^3-) 
B  =  5 — ^  , 


1^  +  162(^-) 
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81  [1-9(4)1 

0      =      2-S-   * 

[l  "  9C$]2 

Hence    the    Integration   of    the   right-hand   3ide   of   equation   (18) 
baoomes    (2) 

i  7(f)2  +   6(f) 

°     360[l-(#($»][x  -$<A>]" 

&3<g)[n*+  162(g)]                   (|) 
J       5 = .    = d(*) 


36o[l-9(V>J  [l-(V)<i>8] 


0 


T 

2    ,2    '    .  „2      .    „.    ~*a' 


a        63  ♦  89KT-)  i  a 


360[l  -  9(V>]         fi  "  ^)(i)2] 


&  M2 

pa      -14U   -   162(^-) 


+  0     360  [1-9(4)  f   [1-3(£>J    d<a> 


4      -    r  _.M2 


S_8i[Ll<oj_  1 

+  0     360[l-9(4)]2'f1-^2     ^ 


.    r    ama  -  27m2    en  r    _  (a!)(A)aj] 

1-240(1   -  I  M2)2  L  ^     a      J 
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L  120(1  -  i  m2)2  m2         S  -  «(*)  Jo 

A 

♦  [      ^*f^2  £n[i-^))]a 
120(1  -  |  irr  ->  a     o 

r     2?(1   -  |  M2)  j_  la 

"^  10(1   -  |  M2)2        [l   -  i(|)]^0 

.      -      U^+27M*)       fa[X.   (^)(f)2J 
240(1   -  |  M2)2  L  b       a 

2 

63  +  891(4-)  75  +  M(f) 

120(1    -  |  M2)2   M2  &  -   M(|) 

+  _m^4[i-i(f)] 

120(1  -  |  i»rr  J  a 

,  27(1  -  I  m2)    r   jjjj     i 

40(1    -   I  M2)2    L    1    -   i(|)J' 
Finally  equation   (18)    becomes 

aRe  240(1   -  |  M2)2  6      a-" 

r^  ♦  M(i) 


63  ♦  89K^-)  p/g  +  M(|) 

+   120(1   -  |  M2)2   M2      n V5  -  M(£) 


2 
,2 


120(1   -  j  M2)2  3  a 
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x 

aR 


27(1  -  \   M2)   r    |(|) 
40(1  -  |  M2)2  4l  -  i(|)] 


2 2 ;^^  !     {i£  +  M2)  ^n  (1  .  X) 

(i  -  |  m2)|         ,        yo"  ♦  M(i) 

1  -  (^)  3   2VE  H    JS   -  M(£) 


♦  (§+  ^)  ^n  [1  -  (\)    %]).  (19) 


Thus  the  velocity  distribution  in  the  entrance  region  for  the 
case  of  uniform  entry  velocity  profile  is  completely  defined  by 
equations  (1)  and  (19). 

When  the  entry  velocity  profile  at  the  channel  entry  plane 
Is  the  one  which  Is  nonmagnetically  fully  developed,  the  same 
equation  (1?)  is  integrated.   Only  In  this  case,  the  boundary 
conditions  are  that  when  x=0,  6=a  and  that  when  x=x,  6=6.   Thus 

1   x       »        7(f)2  +  6<a>  & 

^0^      i  T         J   .  „r  ,  .  o/'l1,     (20) 


|=1360[l-(^-)(f)2][l-i(f)2] 


The  Integration  of  the  right-hand  side  of  equation  (20)  can  be 
performed  in  the  same  manner  as  in  the  case  of  uniform  entry 
velocity  profile. 
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i 

a 


7(£)2  +   6(i) 


1     36o[l-(iV)(|)2][l-i(|)J 


4(4) 


2   uv5 


3va' 

2 


I  3<*fc-)fl4»  ♦  162(^)]  (A) 

;    — IFTS '  *- 


360 fi  -9(^)]  [l  -  (^)(|)2J 


d(A) 


+  / 


a        63  ♦   891(^3-) 
1     360  [l  -  9(!l-)J 


[x.**1] 


4(4) 


M 


♦   / 


a     -144  -  l62(ir-) 
1     360  [l  -  9(^)J2 


+  / 


a      8l[l   -   9<JU] 


^ 
* 


i 1(A) 


1     36o[l-9(^)J2 


V   "  *<*>J 


3va' 


(■     -144 
'  240(1 


•144   -   2?M* 


fe[l"*(i)2] 


2d<i> 


|M2)2   »  L*  "  iT'-i 


63  +  891(4-) 


"'a' 


120(1    -  I  M2)2   M2  LV5   -   M(|)   Jl 

[  ■  aa  t  azag    A  [a  -Wl* 

L  120(1    -  I  M2)2  L  3Va'JJx 


1m2) 


r27(l    -   f  M-)  -,a 

+  I  40(1  -  I  M2)2  '    [l  -  i(£)ri 
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»%:?&■(  fcfr- fry] -*ci-fe) 

120(1   -  J.  mV   M*   *         &  -  K(f)  ,/°"  -  M  ; 

27(1   .  I  m2) 


144: 


♦  -  -2 


« 


tod  -  |  M2)2     [l  -  1(A)]        3 

240(1    -  |  M2)2  M2  J 

2  1    -  T 

2 
63  ♦  891(4-)  ,     [S  +  M(A)](V5  -  H) 


120(1    -   |  M2)2   M2  [S  -   M(A)J  (75  +   M) 

2  .      1    "    (A) 

2" 

3 

27(1   -  \  M2)  |[(|)    -   1] 

40(1    -   |  M2)2         [l    -   1(A)] 


120(1    -  |  M2)2 


2  "  '     I*  "  3*a 
Finally  equation  (20)  becomes 


4(bz 


_x_  144   -   27M2       \  n      1    "  T  (a}     1 

aHe      =      "   240(1    -^M2)2L  ?         J 


)(1    -^M')'L  .  m' 

2  x  -  T 
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or 


.  2 
63  +  891  (t-)       [,/5  +  M(^)l  (V5  -  M) 


120(1  -  |  M2)2  M2     [,/5  -  M(f)]  (75  +  M) 

+    3M  -*■  27M2     ^>n  1  "  baj 

120(1  -  |  M2)2        | 
27(1  -  |  M2)     |  [(f)  -1] 

+  koO.   -  f  M2)2  *  [l  -  i(A)J 


1   -S- 

d~  =  —7~5 — 1 — t-tt  ( -  '¥ +  m2)  ^n  -r3* 

aile     20(2i"T  -  |  -  1  M4)  l    3  | 

3  [(f)  -  l](l  -  \   M2)    11M2  +  ^ 
2Tl  -  ijg)]  2-/E   H 

U-  (f)M](7g-H)     g   M2 

<n  t— ; 1 — ; •»■  (?  +  -Q-)  x 

[V5  +  (f)Mj(V5*  M)     3    2 


i-(^(4) 

fo  a 

1  -  ("g") 


}  •  (21) 


Thus  the  velocity  distribution  for  the  situation  where  the 
velocity  profile  at  the  channel  entry  plane  is  nonmagnetlcally 
fully  developed  is  completely  defined  by  equations  (1)  and  (21). 

In  these  cases,  it  should  be  noted  that  the  nondimensional 
length  parameter  x/afie>  familiar  from  Schlichting's  nonmagnetic 
solution,  is  a  function  of  ft/a  and  the  Hartmann  number  M  only. 
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The  Hartmann  number  can  be  viewed  as  a  measure  of  the  magnetic 
interaction. 

The  Length  of  the  Entrance  Region 
Since  the  boundary  layer  grows  from  zero  until  it  attains 
the  asymptotic  value  and  joins  at  the  channel  center,  when  the 
flow  is  fully  developed  (,a  or  £  =  1.   Therefore,  for  the  case  of 
uniform  entry  velocity  profile,  from  equation  (19)  we  obtain 

f(M  £)|      _  2 (      (16  .  „2.     2 

a  I  |=1     20(2M2  -  |  -  1  «  )  l    3  3 


3 


|(1  -  I   M2)  -  (] 


«2   li. 


y  n  ,  -  (HIT  +  ±?)     ij—  X 

3    2(6)S  M 

(6)s  -  M    3    2  6  / 

Therefore,  the  entrance  length  x  for  the  case  of  uniform  entry 
velocity  profile  Is 


aHe 


a 


20 


-  2   ?2 — rnr  ( -  <¥  ♦  m2)  ^n  I 

-  1(1  -  I  M2)  -  (11M2  +  %    —1     ^n  (6)*  t  « 

3   2(6)*  M     (6)*  -  M 

+  (§  +  if ,  ,n  (1  .  £))  (22) 

According  to  Schiller,  the  entrance  length  is  defined  as  the 
distance  required  for  the  friction  factor  to  come  within  10#  of 


67 

the  final,  fully  developed  value.   From  equation  (1)  we  know  that 
the  velocity  gradient  at  the  wall  will  be 

(2S)     _  u[2  _  2  Z]  2U 

In  equation  (21),  the  minimum  value  of  (£)  is  decided  by  the  fact 

that  the  argument  of  the  logarithmic  function  should  not  be  zero 

or  a  negative  value.   For  instance,  the  minimum  value  of  (£■)  for 
the  Hartmann  number,  M,  of  4  is 

'  a  'iln  =  ~Tt 


»--$-».  (23) 

If  the  velocity  profile  corresponding  to  (£■)    is  assumed  to  be 
the  fully  developed  one,  the  friction  factor  at  the  fully 
developed  condition  for  H=k   is 

<g»  -  fl     =  -IT—'  W 

d*   y=0   at  F.D.  *   'f.D.  (Vo^A)a 

Substituting  equations    (23)    and    (13)    into  equation   (24)    yields 

2u  2k 

(£ii)  m        1   -  JE/12  u        12-75 

ay  y=0   at  F.D.  (Vo"A)a  a  JS 

=     |(       I6       )      =     |(-^~)    •      (25) 
a  12V5  -   6  a  2S  -  1 

The  friction  factor  at  any  other  section  is,  by  virtue  of  equa- 
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tion  (13) 

ay  y=o     6     ad  -  ^f) 

When  the  Hartmann  number,  M,  equals  4,  the  value  of  ^  at  the 
location  of  10^  deviation  from  the  fully  developed  profile  Is 
obtained,  (by  use  of  equations  (25)  and  (26)), 

I  (-£*—)  x  0.9       ^5- 


2VS  -  1  6(1  -if) 


<*>  "  3(|)  -  "  f  '  *#  =  -  1-625  • 


Hence 


.     3  ±  V9  -  6.5 
a  =       2 

=  1.5  -  0.79 

=  0.71  .  (27) 

Finally  substituting  equation  (27)  into  equation  (21),  we  obtain 

xe 

~-     =  0.0185884 

for  the  case  of  the  nonmagnetically  fully  developed  entry  veloc- 
ity profile  at  M=4. 
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The  Pressure  Drop 
Prom  the  Bernoulli  equation,  the  difference  in  pressure  be- 
tween any  two  locations  in  the  entrance  region  Is 

Pl  "P2  =  !U2  "lUl  (2<0 

where  P^^  and  P.,  are  the  pressures  at  xx  and  x2,  respectively,  and 
Ux  and  U2  are  the  velocities  at  the  central  core  at  x,  and  x_, 
respectively. 

From  equation  (13),  we  have 

U,   =  5 U   _      u 

(i  -  -  — )  (i    i  _^i 

11        3    a'  a  ~  J  1T> 

Hence   the  equation   (2^)    becomes 

pi-p2  =  f»2[ hrr h-rJ.      (25) 

The  pressure  drop  at  the  entrance  region  and  at  the  fully 
developed  region  can  be  obtained  as  follows. 

For  x  £  f(M,l),  the  pressure  difference  between  the  inlet 
edge  and  any  location  in  the  entrance  region  is  obtained  from 
equation  (25) ,  as 


r0    r      ~      2 


iu2[ 1 1 1 


(1  .  li)     (1  .  I  fOj' 
3  aJ     U  "  J  T' 
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Since  -* 
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P0  "P=  Z^L   \A.2-  l| 
VJ-  "  3  a; 


If  x  =  x  =  f(M,l),  *  a  1.   Therefore,  the  total  pressure  drop  In 


e   -■.--.-..  a 

the  entrance  region  becomes 


p0  "  P  =  I  a2  •  <26> 
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NOMENCLATURE 

a  Channel  half-height 

B  Magnetic  induction 

E  Electric  field  strength 

J  Electric  current  density 

M  Hartmann  number,  M  =  B  a(§)3 

P  Fluid  pressure 

R  Reynolds  number,  B  =  ka\jf/p< 

u  x-component  of  velocity 

W  Half-width  of  channel 

x.y.z   Space  coordinate 

6  Boundary-layer  thickness 

6*  Displacement  thickness 

f  Fluid  density 

<j  Fluid  electric  conductivity 

8  Momentum  thickness 

//  Fluid  viscosity 

T~  Wall  shear  stress 
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CHAPTEH  IV 

Targ's  Method 
The  principal  difficulty  in  analyzing  the  velocity  distribu- 
tion in  the  entrance  region  Is  due  to  the  nonlinear  inertia  terms 
in  the  equation  of  motion.   Because  of  this  nonlinearlty  of  the 
problem,  all  solutions  obtained  to  date  have  been  approximate 
solutions. 

A  new  method  of  analysis  of  entrance  region  flows,  evolved 
from  ideas  proposed  by  Targ  (Slezkin  (1)),  has  recently  been  de- 
scribed by  Sparrow,  Lin,  and  Lundgren  (Z) .      The  method  involves  a 
linearization  of  the  inertia  terms  by  introducing  a  stretched  co- 
ordinate in  the  direction  of  flow.   A  boundary-layer  model  is  not 
necessary  in  the  analysis,  and  velocity  solutions  are  obtained 
which  are  continuous  in  the  transverse  and  longitudinal  direc- 
tions from  the  channel  entrance  to  the  fully  developed  region. 
This  new  technique  is  applied  by  Snyder  (3)  to  the  investigation 
of  MHD  flow  in  the  entrance  region  of  a  parallel-plate  channel. 

The  geometry  of  the  MHD  parallel-plate  channel  Is  the  same 
as  shown  in  Fig.  1  of  Chapter  I.   A  constant  magnetic  field  B.  is 
applied  In  the  y  direction,  and  the  channel  is  assumed  to  be 
infinite  in  the  x  direction.   The  initial  velocity  profile  Is 
also  assumed  to  be  uniform  and  at  a  great  distance  downstream 
from  the  entrance  the  velocity  profile  Is  the  fully  developed 
Hartmann  profile.  — 

As  can  be  seen  from  Chapter  I,  the  governing  equations  of 
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motion  for  laminar  magnetohydrodynaraic  flow  are 
continuity: 


3u   sv  S 

—  +  —  =  0     or     — 
3x   ay  3x  - 


=  0     or     f-[S   udA  ]  -  0  ,       (1) 


momentum: 

u  <>x  *   v  ay  "    ?  dx  *   ~2  +  H3^  'V^o'-   (2) 

^y 

The  major  assumptions  implicit  In  equations  (1)  and  (2)  are  the 
same  as  those  described  in  the  previous  chapters. 

According  to  the  new  technique,  equation  (2)  is  linearized 
as  follows: 

2 

2     nB    E 

eU)  ufa   =  a(x)  ♦  y^u     _^  c  o  _  u)  (3) 

?y^   "   uo 

where  U  is  the  mean  velocity  defined  as 

U  =  i  /S  udy  (k) 

a  0 

and  £(x)  is  a  yet  undetermined  function  of  x  which  weights  the 
mean  velocity  U,  while  /|  (x)  is  a  second  undetermined  function 
which  includes  the  pressure  gradient  as  well  as  the  residual  of 
inertia  terms  after  the  linearization. 

Integrating  equation  (3)  over  the  upper  half  cross-sectional 
area  of  the  channel.  A,  gives 
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2 
2        rt  B      E 

£(x)U  £-  /  UdA  =   M(x)dA  *  ^J^dAt  -V2  /  (^  -u)dA 

3X  A         A  A3/        °   A   B0 

where,  in  accordance  with  the  mass  conservation  equation  (1) 


ax  /  udA  -  o 


Hence , 


>2„      oB*    En 


0   =  A(x)A  +  V/  *-*  dA  *  -%£  /  (=£  -  u)dA.      (5) 
A  gy^        -T   A  c0 

By  applying  the  divergence  theorem  to  the  integral  Involving 

2    2 
9    u/ay   in  the  above  equation  we  obtain 

»/  ^fdA   =   tfffgd*. 
A  ay  c  *" 

Substituting  this  into  equation  (5)  and  solving  it  for  (x),  we 
obtain 

mx)  -  -f#«ad/*ai!(1.i)  (6) 

c  s'n      ~      UBq 

in  which  the  contour  c  Is  the  circumstance  of  the  channel  half 
area,  au/gn   is  the  normal  derivative  of  the  velocity  evaluated 
along  c,  and  d  is  an  element  of  the  circumference.   For  the 
parallel-plate  channel,  the  only  nonvanlshing  contribution  to 
su/an  is  C8U/9y)y=a- 

A  stretched  axial  coordinate  x*  is  defined  by  the  relation 

dx  =  e  (x)dx*  .  (7) 
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Combining  equations  (3),  (6),  and  (7)  we  obtain  the  linearized 
momentum  equation: 

2  — 
at.    31.  a  0  an  3y±  r  UBq 

or 

2  — 

u^    -    P^-^Sr5       J  +  fsT^(1  --"      (8) 

ax  ay^  a  sy  y=aJ  f  U 

It   is  also  convenient   to   introduce  dimensionless  variables 

He     .     ^  ,  and  M     =     BQ   a(^)*    . 

Hence, 


(9) 


au  u     ,    u   .  u     /aa, 

ax*     =  x*      vfle   a'      =      He   a   'gg' 

9  He   a 


u 

su     =    .f_u  ,g.         u  ga 

ay      "     jl    «  *     »1! 
a 

2  —  —  —       2 

g    U  U   _£_  £&  U      3     3a  U     g   q 

ay2  a  »y  9l  "      a    a  v.  ^      -      fl2  3?2    ' 

Substitution  of   these  gives   the  dimensionless  form  of 
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equation   (8)    as   follows: 

2  — 
Re  a  S^   "   a2L  ^2    V*„1J      J>    U   „' 


or 

Sp   -  %  '  (^7=1  +  w2(1"a)  •  (10) 

The  boundary  conditions  are 

7=1     t     a  =  0  , 

7=0      :    ^     =   0  .  (11) 

3  =  0     :     a  =  1  . 
That  is,  equation  (10)  is  to  be  solved  subject  to  the  no-slip 
boundary  condition  that  u=0  on  the  duct  wall.   In  addition,  the 
velocity  profile  at  the  inlet  will  be  assumed  to  be  uniform 
across  the  section,  I.e. ,  u  =  U. 

Although  the  relationship  between  the  actual  axial  coordi- 
nate x  and  the  stretched  coordinate  x*  is  not  yet  determined,  it 
will  be  convenient  to  set  this  question  aside  until  later  and 
proceed  with  the  solution  of  equation  (10). 

Equation  (10)  is  a  linear  equation  with  constant  co- 
efficients and  is  thus  amenable  to  solution  by  means  of  the 
Laplace  transform.   Defining  the  transform  of  a  as 

a     =  /  ae"*S|3  dg 
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then  we  get 

L(a')  -  sa  -  a(0) 

=  sa  -  1 
in  which 

a(3=0)  =  1  ; 

L(a)  =  a  ;     and 

L(K2)  =  |  M2  . 

Hence,  equation  (10)  is  then  transformed  into  the  form 

or 

«£§  -  (M2  +  8)5  =   (ft     -  ^2  t  3>  .  (12) 

d^r  d7  7=1      s 

To  solve  this  linear  differential  equation  of  second  order,  we 
first  find  by  observation  the  complementary  function 

ac  =  C^W^^-UW   , 

Then  the  general  solution  is  found  to  be 

5  _   „  .(M2*s)*7^  „   -(M2+s)*2 


■fr[#>      -^j.     d3) 

r+s)  L  *"  7=1 


(M"+s)  L  T  7= 
The  boundary  conditions  (11)  can  be  transformed  as  follows: 
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7=1    :    a  =   (/  ae-8**  dg)     =  0      (Ik) 
0        7=1 

-  /*"  2S  e"se  da 

=  o  (15) 

With  these  transformed  boundary  conditions  (\k)    and  (15) .  the 
constants  In  the  solution  (13)  are  to  be  found. 
From  (14)  i  we  have 

c.e"^*  ♦  o-e-^2*8)*  -  -£-[(£)    -  *&■*  1=0. 
1  *■  (mSs)  l  a7  ?=i     s   J 

From  (15)  we  have 

c^M^s)*  -  c2(M2+s)^  =   0 


°1  =   C2  • 
Solving  these  two  simultaneous  equations  for  c,  and  c?,  we  obtain 

_         2 

0  re(M2+s)i  +  e-(M2+s)^  1    [(df^=i  "   s   j 

1  l  J  '         (M2+s) 


r,da,    _  («%s)  ] 

lvd7;7=l  "    s    J 


C,    =   C- 
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Substitution  of  these  constants  into  equation  (13)  gives 

jjdSj    _  (M2+s)jMM2+s)^-M2+s)^  j-.(M2fs)^|C-(M2+3)M 

a  =  _  (M2*s)  [e«^>*  ♦  .-<«2«>*J 

[(g)ycl  -  i^][oosh(MW?  -  ooshuW] 
(M2+s)    [cosh(M2+s)*] 

From  this  we  get 

da     _      [(^^-^JCA.)*^^,)*^ 
d?  (M2+s)    [cosh(M2+s)*] 


d?  ?=1  (M2+s)  [  cosh(M2+s)*J 

Solving  for  (f^)»  ii  we  obtain 

,da>      m     (M2+s)  slnh(M2+s)^ 

**   Y=l  s   [slnh(M2+s)^  -  (K2+s)*  cosh(M2+s)^J 

As  a  result   the   general   solution  of   equation   (12)    becomes 

[ (M2+s)    8lnh(M2+s)*  .  i«f±ai][008h(II2+8)*r 

Ls[slnh(M>s)*  -   (M>3)soosh(M  +s)*3  3        x 


■s[slnh(M Sa)»   -   (M  +s)gcosh( 

(M2+s)    cosh(M2+s)* 

-   cosh(M2->-s)^] 
i[(M2+s)*   sinh(M2+s)*   -   (M2+s)    cosh(M2+s)*] 


[(M2+s)][cosh(M2-s)^  -  cosh(M2+s)^]  (l6. 
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Carrying  out  the  Inverse  transform  C+) ,  gives  a  in  the  following 
two  equivalent  forms: 

1      rr*1°°es&        (M2+s)  rcosh(M2+s)^f-cosh(M2-fs)k]ds 
2TTi   Y-1-  6        s  [(M2+s)5sinh(M2+s)*-(M2+s)oosh(M2+s)*] 

(17) 


or 


.      ^9+1-  2x2eU   -M   \»    (ooshrt-oo8hx)<U 
2TT1    9-1°°  U2-M2)(sinhA.   -  XcoshA.) 


(18) 


in  which 

X   =   (M2+s)*  ,     s   =   \2-M2  .     and    ds  =  2\dA.  . 

The  integrals  of  equations  (1?)  and  (18)  can  be  evaluated  by- 
residue  theorem  (5)  which  states  that  if  f(s)  is  a  single-valued 
function  that  is  analytic  inside  a  simple  closed  path  c  and  on  c, 
except  for  a  finite  number  of  singular  points  a^,  a*,,.,,,   am 
interior  to  the  region  bounded  by  c,  we  have 

m 
/  f(s)ds  =   2TT1     2     Res  f(s) 
c  J=l,s=a 

where  Res  f(s)  denotes  the  residue  of  f(s)  at  a  and  where  the 
integration  is  in  the  counterclockwise  sense  around  c. 

From  equation  (17)  it  is  clear  that  a  simple  pole  occurs  at 
s=0  with  a  residue 
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Hes   f(s)      =      Res  H\ 
s=0  s=0   qv' 

P(0) 
=      qMOl 


I eSB(M2+s)^    [cosh(M2+s)H 

\   slnh(M2t-s)*  -  (M2+s)*oosh(M2+s)i  ♦  § (M2+s)* 

-  cosh(M2t-s)g] _1 


oosh(M Ss)3-(M  +s)~scosh(M  +  s)5-sinh(f 

MtooshM?  -  coshM)  i-,q\ 

-        (sinhM  -  McoshM)  ' 

From  equation  (18),  if  k    is  pure  imaginary,  say  X  =  iy  where  y  is 
real,  it  is  clear  that  poles  will  occur  where 

sinh  k   -  k   cosh  k     =  0 


3inh(ly) 
soshdyJ 


Sl7)t'(   =   iY 
COS1"' 


I.e.. 


tan  Y  .  (20) 


The  roots  of  equation  (20)  give  an  infinite  set  of  eigenvalues  at 
which  poles  occur  in  the  integrand  of  equation  (18).   If  Y   is 
the  nth  root  of  equation  (20),  the  residue  corresponding  to  this 
root  is 

-(Y2+M2)g 
2e  (cos  y„7  -  cos  Y„) 

H_  =  5—5 Ei S_  (2i) 

(Yn+M  )  COS  Yn 


84 


The  solution  for  a  obtained  as  the  sum  of  all  residues  of  the 
integrand  in  the  inversion  integral  becomes 

>shM   -  coshM) 


M(co 
"  "Til 


inhW  -  McoshW 

-(Yn>M2)3 
a,  e  (cos  y    1   -  cos  y„) 

♦  22   2—2 S '  (22) 

n=l        (Yn+M  *  oos  Yn 

The  first  term  of  equation  (22)  is  recognized  as  the  Hartmann 
profile  and  represents  the  fully  developed  solution,  and  the 
series  corresponds  to  the  difference  between  the  actual  velocity 
profile  and  the  fully  developed  profile.   The  limiting  form  of 
equation  (22)  for  M — *0   can  be  written  as 


- ,  T ,    M( coshM f  -  coshM) 

"i,  tt  =  J^L      (sinhM  -  McoshM S 


M — *0 


-lyWlp 


e   n     (cos  Ynf  -  cos  Yn) 


♦  Llm  2M  -      —  j — =— 

M— *0  (y„*M  )  COS  Yn 

in  which,  the  first  term  in  the  right  side  can  be  evaluated 
according  to  L'Hospltal's  rule.   Thus 


M(cog 
(sir 


■^oshMY  -  coshM) 
M_      n      (sinhM   -  McoshM) 


(coshM  1-   coshM)    +   M(?slnhM?-   slnhM) 
„     mQ  coshM   -  MslnhM  -  coshM 

L,         (yslnhMy-slnhM)-)-yslnhM   -slnhM+M(g2ooshMy-ooshM) 
M__Zn  -McoshM  -   sinhM 


T  .    H.    coshM?  -  coshM  +  t    coshMfr  -  coshM  +  I    coshM? 

„   „  -MslnhM  -  coshM  -  coshM 

M 


-coshM  -  M(y3slnhMy-  sinhM) 


Hence, 


■yS* 
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•,     o       »  e     (cos  Y.f-  oos  Y-) 

Llm  a  •  \  (l-f)  +  2    z    o 

M — »>0  n=l        Yn  oos  Yn 

which  Is  identical  to  the  velocity  profile  for  the  case  of  zero 
magnetic  field  dealt  with  in  (2). 
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Stretched  Axial  Coordinate 
(Evaluation  of  £  ) 

The  solution  given  by  equation  (22)  is  in  terms  of  the 
stretched  coordinate  9  and  the  relation  between  9  and  the  physi- 
oal  coordinate  (tfe*)(x/a)  requires  a  knowledge  of  the  relationship 
between  £  and  3.   This  relationship  will  be  found  by  imposing  an 
additional  physical  constraint,  namely,  that  the  local  pressure 
gradient  s p/sx   calculated  from  momentum  considerations  be  the 
same  as  that  calculated  from  mechanical  energy  considerations. 
For  an  exact  velocity  solution,  a  unique  value  of  the  pressure 
gradient  would  be  derived  no  matter  whether  one  employs  momentum 
or  mechanical-energy  considerations.   However,  since  entrance- 
region  analyses  are  necessarily  approximate,  the  local  pressure 
gradients  calculated  on  these  different  bases  need  not  necessar- 
ily be  the  same. 

An  expression  for  the  pressure  gradient  will  be  derived 
first  from  momentum  considerations.   The  momentum  equation  will 
be  Integrated  across  the  channel  half  section.   The  integration 
of  the  Inertia  term  is  facilitated  by  combining  the  inertia  term 
and  the  continuity  equation: 

suu  +  guv  u^u+u2U£U5V 

gx       9y  ax  ax  ?y  sy 

2x       sy'  &x  9f 

where 

9U        av  . 

£7*  i7    "    °  • 
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Hence, 

U2H+V2B     _     »uu       auv  (23) 

ax  sy  9x         ay 

Substituting  this  into  the  momentum  equation  (2)  and  integrating 
the  same  aquation  over  the  upper  half  of  the  channel  gives 

n      a     „  a     Z  a  o  Bn 

=    "  \  !    H  dy  ♦  » '    ^  iy  *  /    -V  (E0  -  uB0,dy 
"  0     9X  0     ay  0       -r 


^  /     u2dy  +   (uv)|  -   (uv)| 

**  0  'y=a  'y=0 


-    "  JS  ♦  »  <T>  ly=a  "   *  ®  ly=0  ^  (Eo  "  UV 

In  which  (uv) [   and  (uv) J   are  zero,  because  u  and  v  are  all  zero 

a  0  , 

at  the  wall  and  v  is  zero  at  y=0;  and  (— )    is  also  zero  because 

»y  '0 

u  is  symmetric  at  the  center  line.   It  is  also  to  be  noted  that 

a     _ 
/  ud  =  Ua.   Hence, 
0 

1*E  „  I  ^/au2dy  -  (£«»,     .  !a!^  (^--D.  (2U) 
f  0X  a  9X  0         a  ay  y=a     /     UB0 

Next,  a  mechanical -energy  equation  may  be  constructed  by  multi- 
plying through  the  momentum  equation  (2)  by  the  velocity  u.   The 
mechanical-energy  equation  may  then  be  integrated  over  the  upper 
half  of  the  channel. 
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a     i  a 

J     u     -dytj     uv  -  dy 


a     n_  a     S2V  a  a„Bn 

2 

ay 

In  which 


t      a  a     ,2  a  a  B. 

?  /  u  |f  dy  ♦    *J   u  2-4  dy  ♦   J     -a-2  (E  -uB   )udy 

f  0       9X  0      ay  0       J  u       u 


jVgdy    =     ;a|^dy./au^dy 
0  "x  0      s  0  * 


3     ,         3.  ."        S\X~ 

=     -  J      u^dy  -  /     u  -sv  dy 


a        ^2 
3X 


r  3U     ,  f      ?U  V    j  r  SUV     . 

;     vu  -^  dy     =     /     -ff-  dy  -   /     u  — -  dy 
q  93  0  y  0  y 


-i    i  a  a 

2  p  <9  uv   „ 

=    u  v      -  ;    u  T7  dy 


J   ■»•» 


Hence, 


r        2  8U   J  ,  au    . 

/     u     rr  dy  +   /     uv  -^  dy 

=     £  /S  u^dy  -  /"  u  f£  dy  .  /*  u  222  dy 


0  0 

_      a      .,  a  2 

3     f         3j  t         /9U  suv>  . 

-   a~  '    u  dy  "  '    u(^r  +  77)dy 


—  /     u3dy  -   /     u(u-+v  — )dy 
pi  J  '        q        v      ax  ay'    ' 


?y 
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&      p      XL      j 


.  I  ;a  u  ££  dy     =     .  i  ££.  /a  udy 

because    (u)  =   0   and    (|^)y=0   =   0. 

2 
a  J  fl.  a  B.E.      a  aB.      a     „ 

/     -VA  (E     -  uB0)udy     =     S-2S.  J     udy  .  -A  J     u2dy 
0-Tuu  JO  r0 

oBnEn      _  OBn      a      . 

=     -P  (Ua)    -  -£  I      u2dy 
J  JO 

CTB2   U      EAa 


0   "    ,"0"        lf"     2.    . 

°        Do       u  o 

Substituting  this  back  we  may  obtain  the  pressure  gradient  from 
the  mechanical  energy  equation  as 

.a  „3      #    .a  „„  2 


"pJx  =  w  Sx  *     -5-  dy  ♦  —  /  (§-)  dy 
J  ,9X     Ua  3X   0   ^       Ua  0   ay 
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aB2   U       E, 


£—  (r5 h?  I    u2dy)   •  (25) 

"  UBQ        alT   0 

The   pressure   gradients   given  by  equations    (2^)    and   (25)    may  now 
be   equated   to   give 

I  ±    rS  „2.  ,i>3u,  °eB0   "    ,  E0         ., 

2 

Ua  a      0  Ua  0  * 

aB0   I   ,  E0  1      ;a     2,    . 

5~  (" T?   J      u  dy)    . 

J  UBQ        aU     0 

Replacing   (d/dx)    by   (— )(d/dx*)    and   rearranging,    we   get 
e 

in*  <£***-$($«> 

e  = *-§-— 3-5 (26) 

,->  n       a  2  a  B„   U  a   B„   U     a     2 

/*     ui  "     f    iau\    a„  e   0  e   0         f      u      . 

(r  ~7)  v-a  *   ~   J    (»vJ    dy P +       Fa —   J      ^2   dy 

y  y  a        Ua  0        y  J  fa0U 

in  which,    from  equation   (9),    we   have 

_£_  9  1  _1 3_ 

<9  x*     =  x*     "    B  a     =      H  a  a  0    • 

Ha  e  e 

e 

U2     =     a2   U2    . 

U3     =     a3  P   . 


dy     =     ad  7 
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u 

0  ~  —               — 

£u    u  u    u  ffa 

sy  3  j  I  a  =  a  a^  • 


and 


aeB0  U     JeB0a2   UP    M2  ff )> 
T  =    ^    '  a2        a2  * 

Substituting  these  into  equation  (26)  we  obtain  the  expression 
for  6   in  terms  of  dlmensionless  variables. 


1 L    "•1     2-2     ,,  „      1    fX  a3  U3 

2  <?x« 


(/     a2   U2  ad?  -  1  /     SLjL.  ad^) 


B  a     "A        0  U  0 

f(M,3)  = r— 2 -» . = 

(r-S^l        +  —  '  -p  •  a  /  (jzO  d7*-2  M     /  (a  -i)d7 
a     a  <?7  7=1      Ua     aT  0*7  eT  0 

= r5 — g j .  (27) 

(f#>  *  /   <57>   d?  +  M2  /   (a2-l)d>Z 

l  7=1      o    '"  o 

The  right  side  of  equation  (27)  is  a  known  function  of  9  from 
equation  (22),  and  thus  the  variation  of  £  with  @  is  specified. 
It  is  clear  that  the  (£-3)  relationship  will  also  depend  on  the 
Hartmann  number  M.   The  physical  Coordinate  x  is  related  to  the 
dimensionless  stretched  coordinate  3  by  the  following  inte- 
gration: 

From  equation  (7), 

dx  =  £■ dx* 
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where 


Hence, 


dx*  =   Ha  dg  . 


X  0 

/  dx  =  Ha  /  e dg 
0  e   0 


8 
Ha  /  6   dp  .  (28) 


Equation  (28)  can  be  evaluated  numerically. 

Equations  (22),  (27),  and  (28)  thus  give  the  solution  for 
the  velocity  profile  in  the  entrance  region.  From  this  basic 
solution,  the  calculation  of  other  physical  quantities  of  inter- 
est is  possible. 
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Pressure  Drop 
The  pressure  distribution  may  be  obtained  by  integrating 
either  the  momentum  equation  (2*0  or  the  mechanical-energy 
equation  (25).   The  momentum  equation  has  a  somewhat  simpler  form 
and  will  be  employed  in  the  following  development.   The  pressure 
drop  between  the  duct  inlet  (pressure  P_)  and  any  axial  location 
(pressure  P)  may  be  determined  by  Integrating  equation  (24)  from 
x=0  to  x=x. 


2  7T 


E 


0 


,i  x  a(u-u-)  x  a  Bn  U 

-  /   (- *—)  dx  +  /  S-§—   (-*-  -   l)dx 

a0      "    y=a      0     /     UBQ 


P0_P       ■*'    1       *       -a   '   1       >>  -X  .»"> 


a  i    s1  y   y=a 


0     "  '    '  —         J     UBQ 

where  the  identity 

u  =  u.  +  (u-u.) 

has  been  used  in  evaluating  (au/ay)    .   Employing  the  dimension- 
less  variables  expressed  by  equations  (9)  we  have 


/X     i  {=  ^  u2dy  dx 
0  ax    0 


1  /     aU2  a2  d7 


U2  /     a2  d  7  , 


i>   ,x  3Uf 

Ij<77>y=ad* 


=    J»  (_JL)        S  .  x 
a     0  i  y=a  a 


a2       *7    ?=1 

i72     tc     3af 


x    a(u-\xr) 


r  /     (- 


)  ,  .   dx 


p  y      y=a 


|/(-JL-JL-l         «     ad9 


3(f)' 


y=a        e 


o  B       3(a-a   ) 


9^ 


and 
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f    UBQ  a'   UBQ 

He  a  UBQ 

Substituting  these  into  equation  (29) .  we  can  obtain  the  pressure 
distribution  in  dimenslonless  form: 

pq  : F  _  i.  (*)/  M2(  _  ^o_}  _  (^f}    \  + 


f  i  o  e       (o-ou)  > 

+  2(/a8*?-  J  6  i-sf-)tml   dp) 


(30) 


where  the  first  bracket  term  represents  the  pressure  drop  that 
would  result  from  a  fully  developed  flow  and  the  second  bracket 
term  represents  the  correction  due  to  the  entrance  flow. 
Since,  from  equation  (22), 

MfcoshM?  -  coshM) 
af  "   (sinhM  -  McoshM) 

("?3")„  i   can  t>e  evaluated  as 

l90t\  M2     sinhM 

v«?f  'i=l     ~  McoshM  -  sinhM  ' 

And  by  defining 

UB 


1    =  r^ 
"Bo 

and 
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It       {    1     3       3   •?(a-af.)    , 

Ku.ifj  =  2(;  .v-^iy-'tJ1'  • 

Equation  (30)  can  be  written  as 

£s_l4  _  Mi(*jf(i.rf)4 jiaia _) 

1  (>jj2      He  ka;\  *x  v>    *   MooshM  -  slnhM  / 

The  coefficient  of  (-^-)  (J )  In  equation  (31)  Is  seen  to  be  a 
Ke  a 
function  of  M  and  0   and  may  be  defined  as  the  fully  developed 

friction  factor.   Thus  we  may  write 

IT^       f  t  (i>  *  K(w'  t  •*  (29) 


k?v 


where 


2 


f      =     2M' 


{(1-0) 


slnhM 


McoshM  -  slnhM 


From  this  it  can  be  observed  that  it  is  the  fully  developed 
friction  factor,  f,  not  K,  that  is  influenced  by  the  electric 
field.   The  term  K  may  be  considered  as  a  correction  term  to  be 
added  to  the  fully  developed  pressure  drop  to  account  for  the 
entrance  effects. 
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NOMENCLATURE 


Hartmann  number,  M  =  B-  a(^)* 


A  One-half  the  channel  cross-sectional  area 

a  Channel  half -height 

B  Magnetic  field  intensity 

E  Electric  field  strength 

K  Pressure  correction  term 

M 

P  Fluid  pressure 

H  Reynold's  number,  R  =  fUa/^ 

s  Laplace  transform  variable 

u  x-component  of  velocity 

v  y-component  of  velocity 

U  Mean  velocity 

x,y,z  Space  coordinates 

x*  Stretched  axial  coordinate 

a  3  =  x»/Rea 

v  Eigenvalues  of  tan  v  =  v 
n  n   'n 

£(x)  Scalar  factor 

7  Dimensionless  transverse  coordinate,  7=  * 

'   a 

a  Dimensionless  x-component  of  velocity,  a  =  — 

U 

aQ  0-   evaluated  at  y=0 

a.  Laplace  transform  of  a 


A  A,  =  (M2+s)* 

A  (i)  Defined  by  equation  (3) 

M  Viscosity 

V  Kinematic  viscosity 
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Density 


o       Electrical  conductivity 


0   =  E0/UB0 
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Recently  attention  has  been  directed  toward  the  effect  of  a 
magnetic  field  on  the  flow  of  an  electrically  conductive  fluid, 
usually  referred  to  as  magnetohydrodynamics  (MHD) .   MHD  is  im- 
portant because  it  finds  many  applications  in  engineering,  such 
as  an  MHD  generator,  the  pumping  of  conducting  fluids,  plasma 
confinement  for  the  fusion  reaction,  propulsion  and  flight  con- 
trol for  rocket  and  hypersonic  aerodynamic  vehicles. 

In  many  MHD  applications,  the  flow  of  the  extremely  high 
temperature  fluids  is  seldom  fully  developed.   For  this  reason,  a 
study  of  the  velocity  fields,  boundary  layer  development,  and 
friction  factors  in  the  entrance  region  of  MHD  channel  is  of 
practical  importance  and  has  been  a  subject  of  investigation  in 
recent  years.   Although  the  exact  analytical  solution  is  not 
available,  many  approximate  solutions  have  been  presented.   The 
approximation  methods  of  solution  can  be  classified  into  four 
different  categories:   the  momentum  integral  method,  the  linear- 
ization method,  the  matching  method,  and  the  finite  difference 
method. 

The  purpose  of  this  report  is  to  study  laminar  MHD  flow  in 
the  entrance  region  of  a  flat  duct  by  presenting  in  detail  the 
solutions  according  to  the  momentum  integral  method,  the  linear- 
ization method,  and  the  matching  method. 

A  brief  introduction  of  basic  governing  equations  of  MHD 
flow  —  Maxwell's  equations,  the  continuity  equation,  and  the 
modified  Navier-Stokes  equation  —  is  given  in  Chapter  1. 


In  Chapter  2,  Schlichting ' s  matching  method  applied  by  Roidt 
and  Cess  to  MHD  entrance  region  flow  is  presented  in  detail.   The 
flow  field  is  divided  into  two  sections  and  an  appropriate  analy- 
sis utilized  in  each.   In  the  section  near  the  inlet  a  boundary- 
layer  formulation  of  the  equation  is  used  and  a  solution  developed 
in  a  series  stream  function  with  Blasius'  function  as  coefficients. 
When  this  solution  becomes  unwieldy,  an  exponential  velocity 
deviation  from  the  fully  developed  flow  is  assumed  and  this  is 
joined  to  the  boundary-layer  solution  to  complete  the  descrip- 
tion of  the  flow. 

Chapter  3  contains  a  detailed  account  of  Schiller's  integral 
method  applied  by  Maciulaitis  and  Loeffler  to  the  problem.   A 
momentum  integral  equation  is  derived  first.   Then  by  applying 
the  assumed  velocity  profile,  which  involves  the  boundary  layer 
thickness,  to  the  momentum  integral  equation,  a  relation  between 
the  boundary  layer  thickness  is  obtained.   The  solution  is  com- 
pleted after  carrying  out  the  integration.   Two  cases,  the  de- 
veloping of  velocity  from  a  uniform  profile  at  the  entry  to 
Hartmann  velocity  and  the  developing  of  velocity  from  a  parabolic 
one  at  the  entry  to  the  Hartmann  one,  are  considered. 

In  Chapter  4,  the  linearization  method  evolved  from  ideas 
proposed  by  Targ  and  applied  by  Snyder  is  presented  in  detail. 
The  nonlinear  partial  differential  equation  of  motion  is  linear- 
ized by  the  assumption  that  the  convective  terms  are  a  function 
only  of  the  direction  of  flow.   A  stretched  axial  coordinate  is 
introduced  in  the  linearization  of  the  convective  terms  is  the 
key  point  of  the  method. 


